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double layer metal, and (4) the investigation of 2-D modeling of MOSFETs. 

Using the 2-D diffusion program developed in the previous contract period, Simula** 
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missing equipment being the SiO ? vapox reactor. I 

The double layer metal activity initially emphasized wet chemistry techniques. By [ 
incorporating the following techniques: (1) ultrasonic etching of the vias, (2) premetal 
clean using a modified buffered HF, (3) phosphorus doped vapox, and (4) extended sinter- 
ing, yields of 98% were obtained using the standard test pattern. 

A two dimensional modeling program has been written for the simulation of short 
channel MOSFETs with nonuniform substrate doping. A key simplifying assumption used is 
that majority carriers can be represented by a sheet charge at the silicon dioxide- 
silicon interface. The program is not complete. In solving current continuity equation 
the program does not converge. However, solving the 2-dimensional Poisson equation for 
the potential distribution has been achieved. The status of other 2-D MOSFET simulation 
programs are summarized. 
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SUMMARY 

During the reporting period of this contract, four activities 
were in progress: 

* 

1) the completion of the 2-D diffusion studies in SOS, 

2) the set-up of a sputtering system, furnaces, and 
photolithography related equipment at MSU, 

3) experiments on double layer metal, and 

4) the investigation of 2-D modeling of MOSFETs. 

Using the 2-D diffusion program developed in the previous 
contract period, simulations of various phosphorus and boron 
diffusions in SOS have been completed. 

The MSU I. C. Laboratory has in place the furnaces, the 
photolithography facilities and metallization system to reproduce 
the NASA/Huntsville double layer metal process - the only missing 
equipment being the Si0 2 vapox reactor* 

The double layer metal activity initially emphasized wet 
chemistry techniques, By incorporating the following techniques: 



(1) ultrasonic etching of the vise, (2) premetal dean ualng a 
modified buffered HF, (3) phosphorus doped vapox, and (4) extended 
sintering, yields of 982 were obtained using the standard test 
pattern, 

A two dimensional modeling program has been written for the 
simulation of short channel MOSFETs with nonuniform substrate 
doping. A key simplifying assumption used Is that the majority 
carriers can be represented by a sheet charge at the silicon 
dioxide-silicon Interface. The program is not complete, In 
solving current continuity equation, the program does not converge 
However, solving the 2-dimenslonal Poisson equation for the 
potential distribution has been achieved. The status of other 

r 

2-D MOSFET simulation programs are summarized. 



t; REDISTRIBUTION DIFFUSIONS FOR ION-IMPLANTED 
PREDEPOSITS OF BORON AND PHOSPHORUS IN SOS FILMS 

The objective of this work was to produce curve* describing 
the variation with diffusion tine and temperature of the junction 
depth, sheet resistance and integrated inpurl ty dose. This data 
has been generated for boron and phosphorus redierributed in 
nitrogen, dry oxygen and steam ambients for <111> oriented SOS 
films. The following section presents discussions of the implan- 
tation and redistribution model, further program develop, the 
computational procedure and of the computed results. 

A. The Redistribution Model 

There are three aspects of the redistribution model which are 
considered: (a) the implanted profile, (b) the oxidation model, 

and (c) the dlffusivity model. 

a.) The Implanted profile, The implanted profile is described 
by the Gaussian function, 

(y - R > 2 

°(y) ■ aR P ) < 1 > 

P 

where C is the concentration, y is the distance from the entrant 

silicon surface, R is the range and AR^ is the straggle for 

P ' P 

the implant. The peak concentration, c max » is related to the 
implant dose, Q imp by: 


Redistribution data has been generated for the following condi- 
tions: [1,1] 


* 

imp ‘ 

5xl0 12 , 10 1 ’ 

1 _ , „13 .^14 -2 

, 5x10 p 10 cm 

v 

0.2735 pm 

80 ksV boron implant 

P 

0.1727 pm 

150 ksV phosphorus. 

AR s 

0.0665 pm 

80 keV boron Implant 

P 

0.0440 pm 

150 keV phosphorus. 


The dose* are light to noderate resulting in concentrations no 
X8 

heavier than 6xlO xo cm , and the range-straggle values are typical 
of those employed at MSFC. It is assumed that all of the ions 
become activated shortly after redistribution begins and thereby 
diffuse by a substitutional mechanism Involving vacancies, 

b.) Oxidation model; The oxidation model is assumed to be the 
Same as for bulk silicon and the data of Deal et. al, [1.2] has 
been used to calculate the oxidation rate according to: 

dx 

* B/(2x 0 + B/C) , ( 3 ) 

where B and C follow Boltzmann-like temperature dependences. 

Figures ( 1 ) and ( 2 ) illustrate the oxide thickness dependence 
upon time and temperature for both dry Oj and steam ambients. 

During the oxidation, the silicon film thickness is reduced 
according to: 

(4 ) 


W 


- W - ax„ 
o o 


where W 0 Is the Initial file thickness, taken to be 1 pm, and 
<* « 0.45 la the ratio of the densities of S10 2 to silicon# 

Redistribution date is given for W Q m 1 ye and an initial oxide 

; • 

thickness of x_ m 300 A. 
o 

c.) Dlffuslvlty model: The dlffuslvlty model for boron was 

discussed in an earlier report [1.3] and it includes a linear 

dependence of the dlffuslvlty upon the vacancy concentration as 

well as the field-enhancement effect* The dlffuslvlty model for 

phosphorous Includes only the field-enhancement effect which is 

sufficient to describe the non-linear behavior of phosphorous 

19 -3 

diffusions at concentrations lower than 10 cm as shown by 
Barry [1.4] aMd Fair and Tsai, [1.5] . The dlffuslvlty temperature 
dependence is after Fair [1.6] and Fair and Tsai [1.5] adaptation 
of data by Ghostagore [1.6] . For either boron or phosphorous the 
effective diffusivlty is given by: 


D eif * D(u) u / S u 2 + ! >» ( 5 ) 


where, 

u “ C/2 n^ , ( 6 ) 

and* 

D(u) * D* u , for boron, (7 ) 

“ D* , for phosphorus* ( 8 ) 


and where n^ is the intrinsic carrier concentration at the 

diffusion temperature and D* and D* are the intrinsic diffusi- 

s - - , P 

vitles of boron and phosphorus: 


D* • 3.17 «p (-3.59«V/ k,T) a’/nc., 

D* - 3*83 «p (-3.66.V/ kjT) ( 9 ) 

B. Further Program Development 

The program which was used to generate, the data haa been 
described in detail in an earlier report. It was noted that the 
program was developed in such a way that one could take advantage 
of a normalization procedure for predeposition diffusions end 
generate data applicable to different film thicknesses. However, 
ibis not possible to gain such an advantage for redistribution 
diffusions involving ion-implants or growth of an oxide. Then 
the program was used to generate data, it was discovered that some 
other features of the program are extraneous unless further refined 
The program was developed to account for both thin and thick 
oxides such as would be encountered in some practical situations. 
However, such a simulation requires the incorporation of a warped 
grid system, a modification which would require considerably more 
effort. Therefore, the variable oxide feature is of limited value 
at this time, since the program, at best, only approximates the 
conditions for growth of a very thin oxide during redistribution. 

A modification was made which allows accurate treatment of 
redistribution under oxidizing conditions when only a single oxide 
thickness is Involved. The original program treated the oxidation 
process with regard to the boundary conditions; however, unlike the 
case of bulk silicon* one must also account for the reduction of 


4 



the silicon film thickness. This feature is now included in the 
program* During the simulation of a redistribution in an oxidising 
sab lent, the vertical grid spacing continuously shrinks while the 
horizontal grid spacing is constant. The aodificetlon does not 
■how up on logic flow diagram* at the level of detail which has 
previously been given* For completeness, s new lisfciag of the 
effected main and sub-programs in given in Appendix C. 

C. Computational Procedure 

The progrea described in «n earlier report, end modified as 
described in the preceding, wee used to generate the date. Two- 


dimensional data was obtained in the form of Isoconcentratlon con- 
tours for typical situations. The bulk of the data which can be 

* 

correlated with experimental measurements is generated using a 
quasi-one-dimens ional model in a manner described in a previous 
report. [1.3] A brief review of the procedure la given in the 
following. 

For generation of sheet resistance, junction depth and inte- 
grated Impurity dose data as a function of time and temperature, 
only a one-dimensional profile need be calculated. This is accom- 
plished by making the horizontal grid only three units wide but 
keeping the field six film thicknesses wide. Periodic boundary 
conditions for the horizontal dimensions are employed in the 


program and result in a calculation which produces the vertical 
profile equivalent to a none-dimensional model. Thus without 
sacrificing the generality of the program for treating 
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two-dimensional cases, the amount of computing time Is drastically 
reduced when the data that Is desired does not require the full 
power of the program. 

The vertical grid varies frpm thirty one to sixty one points 
as required for accuracy in details of the profile, and most of 
the data is not sensitive to the number of grid points used if the 
number is chosen in this range. For the purpose of illustrating 
the unusual nature of phosphorous profiles, the larger number of 
points was required. 

D. Discussion of Results 

First, some of the unusual behavior of redistribution diffusions 
in SOS films will be discussed in this section. Next, the format 
for the calculated curves will be discussed, and, finally, the 
bulk of the generated data is given in Appendix A and B without 
further comment. 

Figures < 1 ) and ( 2 ) illustrate the oxide thickness growth 
and silicon film thickness reduction as functions of time for 
<111> silicon films oxidized in steam and dry 0^ ambients, the 

■ o 

evolution traced beginning with an initial oxide of 300 A thickness 
on an SOS film of 1 pm initial thickness. The curves are shown for 
four temperatures. The data are necessary for interpreting some 
of the results for simulated redistributions. 

Figures <3 ) and ( A ) show impurity profiles for boron and 
phosphorus implants being redistributed in a steam ambient at 
1000 deg. C. The profiles are all plotted with a common origin as 
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would be the case for experimentally derived profilea where the 
Sl-SK^ interface would serve aa the logical origin. However , 
the profilea are normalized with respect to the film thickness 
which is of course shrinking. The boron profiles are not unusual 
but show the wellknown leaching effect due to segregation into the 
oxide. The phosphorus profiles show the effect of impurities 
being rejected from the oxide. There is a pile-up of impurities 
in front of the advancing Si-Sit^ interface and then a dip which 
eventually disappears. It is easy for one to draw an erroneous 
conclusion from observing the profiles, because it appears that 

the integrated dose should increase or at least remain constant 

»* 

and the sheet resistance should decrease with time. This is not 
true. Although the segregation coefficient favors phosphorus in 
Silicon vs. SiOg, eventually all of the phosphorus will be in the 
Si0 2 when the SOS film is completely oxidized since the model 
assumes that there is no diffusion into the sapphire. 

Figures (5-7) Illustrate the behavior of the junction 
migration, sheet resistance variation, and integrated impurity 
dose variation over a long period of time. All of the curves are 
plotted with respect to normalized time, and true time is obtained 
by multiplying by the normalizing time value given on the plot. 

Junction depths are in microns, sheet resistance values are in ohms 

_2 •' /' .. ■ 

and dose values are in cm units unless otherwise marked. The 

curves are given in the typical format for all of the data. 

For an ion-implanted profile, there are in fact two junctions 
until one of the junctions emerges at the Si-SiO^ interface. 





TEHP * 1000.0 

NORMALIZING TIME -114. 1 HIM 
AMBIENT * STEAM 
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Therefore, the sheet resistance values ere for the burled layer 
until the front Junction disappears* This typically occurs in a 
short tine compared with that for through-diffusion of the back 
junction* Figure ( 5 ) Illustrates the through-diffusion of the 

back Junction for the heavier doses. This always occurs for 

* 

redistribution in a nitrogen ambient but not necessarily so for an 
oxidizing ambient. After the through diffusion, or even before 
for light doses, the junction depth will eventually decrease due 
to the reduction of the film thickness or due to the relatively 
slow advancement of the junction with respect to the moving 
Si-Si (>2 interface. In some of the data presented in the appendix, 
the junction appears to remain almost stationary for this same 
reason. The variation of the sheet resistance and dose with 
redistribution time also may appear strange when compared with 
results for bulk silicon; however, consideration of the previously 
mentioned factors also explains these results. 

Two-dimensional isoconcentration contours are given in the 
appendix for the various ambients and the two Impurity types. The 
results are not as remarkable as those given in the last report 
which were for chemically predeposited boron. In that case, 
there was initially a heavy concentration of fast diffusing 
Impurities at the Si-SiO^ interface which were strongly retarded 
due to the segregation phenomena. This does not happen with the 
ion-implanted predeposit because the initial profile lies below the 
interface at which the segregation phenomena is effective. 
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II. DOUBLE-LEVEL METALLIZATION TECHNIQUES 


At Long Tern Objective 

Integrated circuit; technologies which incorporate two -level* 

5V 

of metallization into the interconnection scheme are often limited 
in density and yield because of limitations in the double-layer 
metal technology itself. With a dual level aluminum and deposited 
(low temperature) silicon oxide as the base system, the primary 

... .jj . 

objoctive of this program is to evaluate the basic yield limitations 
of that system and Investigate alternate approaches in order to 
enhance the yield and/or facilitate Improved circuit density. With 
the pressure of device (or design rule) scaling as an on-going 
philosophy, high yield for a large number of small visa and small 
lines becomes extremely important. 

B. Fabrication Facilities at 
Mississippi State University 

Most of the facilities work during this contract period was 
directed toward the development of a sputtering system for pre- 
paring aluminum and aluminum-alloy films in support of the double 
layer metal program. A photograph of the completed system is 
shown in Figure 3. Briefly, the system consists of the following. 

The basic vacuum system is a Varian-NRC 6 inch oil diffusion 
and mechanical roughing pump equipped for automatic and manual 
operation. The diffusion pump is equipped for LN^ cooling of the 

17 
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Figure 8. Photograph of metallization system 
using sputter-gun source. 
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cold trap. Tha sputter gun and power supply vara obtained 
iron Sloan Tachnology. Tha •piittarlng chamber van designed and 
built at Mississippi Stats. It consists of a Corning 12x18 glass 
cylinder and an alunlnun top plate machined to accomodate up to 
three sputter guns. A cold cathode discharge gauge was constructed 
and Installed In tha baseplate of the chamber to Measure the 
pressure during the sputtering operation. A throttle value with 
several threaded holes for accomodating plugs is operated by one 
of two mechanical feedthroughs in the baseplate. A lift mechanism 
vith a reversible motor was designed and constructed for raising 
the top plate and sputter guns. 

Inside the chamber Is equipped with a rotating table which 
accomodates up to eight wafers of l%-2" In diameter. The table is 
driven by a vacuum sealed shaded pole motor through a magnetic 
coupling at 7 rpm. The entire motor-table assembly Is rotated by 
a chain-sprocket-mechanical feedthrough arrangement through three 
sputter-gun and two mask positions! The mask is attached to the 
rotating assembly and provides one hole through which the sputter- 
gun deposits metal on the wafers. A crystal film thickness sensor 
is located beneath the sputter-gun and receives a deposit through 
a hole at an unused wafer position on the table. 

The sputtering chamber opens into a class 100 clean bench in 
order to maintain a high level of cleanliness. The system is 

located in the metallizing and bonding room of the microelectronics 

* * * 

•* 

laboratories in Slmrall Engineering Building, and this room was 
designed with air conditioning and filtering units to maintain a 
class 10,000 environment . 


The Installation of a six-tuba Thermco Ranger diffusion 
furnaca was completed with tha addition of a venting system for 
exhausting tha acavangar boxas. 

Complete PMOS circuits hava bash fabricated In tha MSU 
facilitlaa. 

A summary of the Microelectronics facility at HSU is provided 
in Appendix E, 


C. Double Layer Metal Experiments 


Reviewing the literature has disclosed many investigations 
relative to the surface oxidation processes of aluminum, including 
recent reported results by Bell Lab personnel which directly * 


relates to double layer metal processing. These studies show Chef 
the aluminum surface is extremely reactive and grows Al^O-.nH^O 


films rapidly in water based solutions to thicknesses which 
cannot be reduced during the sintering operation. Even the growth 
rate of A^O^ on A1 in air at room temperature is reported to be 
in the range of 30 \ , This compares with the Bell Lab’s estimate 
of 55 A of AlgO^ which can marginally be reduced by sintering. 

After careful review of the basic double layer metal proceSB, 
certain ones stand out as suspect: (1) The via etch using buffered 

etch, (2) the Initial oxidation during vapox deposition, and 
(3) the photoresist removal. 

Consequently, three sets of experiments were initiated: 


1) modification of the wet chemistry techniques, 
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2) dry chemistry processing end/or new insulation materials, 
e.g. , polymide, nitride, 

3) the incorporation of in situ hack sputtering » 

Dramatic Improvements in yield vers achieved vlth relatively 
minor changes in the vet chemistry processing. The effective yield 
in the number of vias having a total contact resistance of less 
than 500 ohms has been increased to the 98-99 per cent level for 
a test pattern having 560 vias of varying sizes by the use of 
vet chemistry processing techniques and post-heat treatments. 
Alternations in the process vhich accounts for these results 
consists of the folloving: 

1) etching the vias (intermetal contact windows) through 

* 

C.V.D. - SiOg using an ultrasonic etch process in the dark has 
rendered improved results over that vhereby etching vas accom- . 
plished by dipping the vafer in stirred B.O.E. 

2) chemically cleaning the first level metal through the 
vias vith 1:1:1 ethylene glycol: buffered HF^O solution for 
20-30 seconds just prior to depositing the second level metal. 

3) the use of phosphorus doped vapox to allow for faster 
via etch times and the ability to perform post heat treatments 

of temperatures in the 500 °C. range without the cracking (crazing) 
observed vith undoped vapox. The use of phosphosilicate glass also 
tends to give a lover initial contact resistance .over that obtained 
vith undoped vapox. 

4) the use of post heat treatment (sintering) in the vicinity 
of 500 c C. tends to drastically decrease the total via contact 

21 


resistance thereby increasing yield* On wafers exhibiting initial 
poor yields, sintering can improve this yield by 500-800 percent. 

On wafers having good yield initially, the contact resistance can 
be substantially improved. (Average contact resistance of 150 and 
200 ohms have been obtained for 560 vias of 0.5 and 0.2 mil square 
vlas , respectively . ) 

The details of these double layer metal experiments relative 
to wet chemistry have been reported in a special report, but are 
included for convenience in Appendix F , "Sintering Effects on 
Double Layer Metal for VLSI Applications". 

The exploratory processing involving dry chemistry required 
the use of plasma etching and plasma deposition equipment not, 
available at either the Marshall Space Flight Center or Mississippi 
State University. Consequently, the logistics involved in pro- 
cessing the wafers have been difficult. 

Four types of insulators were to be investigated: 

. ■*' 1) 5 % phosphorus doped silox, 

2) plasma nitride, 

3) the GAF version of polyimide, and 

4) the Hitachi piq. 

Whenever possible, dry chemistry processing was to be compared with 
"standard" wet chemistry. The wafers were processed through the 
first metal mask at Marshall Space Flight Center in standard 
fashion. The polyimide material was processed at American Micro- 
systems, Inc. (AMI) Using the standard recommended procedures by 




the manuf acturees , The nitride and sllox were deposited at AMI 
using their standard processes, the former on an Applied Materials 
plasma reactor. The wet chemistry processing was performed at 
AMI using AMI'S standard processing. The plasma etching was 
performed by Davis and Milder on one of their production parallel 
plate plasma system using their proprietary processes. 

The wet chemistry split involving sllox was lost when the 
standard "pad" etch attacked the metal completely. The Same 
batch etchant continued to provide good results on AMI material, 
Since the first layer metal at NASA was deposited via sputtering 
and with many different impurities, compared with approximately 
1% Si:.Al deposited with electron beam, the difference in the 
metal film is believed to be the source of the problem. 

The plasma nitride material and sllox material, deposited by 
AMI and plasma etched by Davis & Wilder, was processed Without 
any problems other than long delays. This material plus the 
polyitnlde material has finally been, received for the final 
processing here. Although, no electrical results are available 
at the present, the microscopic examination displays the expected 
sharper edges of the dry processed material vs. the standard wet 
chemistry. 

In the meantime, the results of the in situ back sputtering 
experiment were obtained. As was expected, the in situ back 
sputtering split gave much improved yields over that for the 
standard material - the detail results being published elsewhere. 

The results of these experiments points to the requirement to 
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very effectively remove the Ai,^ prior to the second metallzation 

or minimize its formation during wafer processing to within the 
6 

order of 30 A which a good sintering treatment can reduce. The 
in situ cleaning forgives prior processing which tends to grow 
the AlgQj such as photoresist removal and buffered HF etching 
of oxide. 

The dry processing has the potential in the future of 
resolving smaller via sizes and, simultaneously, avoids the 
harmful H^O to accelerate the aluminum oxide (and hydrate) 
formation. Although dry processing to remove the photoresist is 
routinely done one should expect AlgO^ formation due to the 
presence .of active oxygen in the "ashing" process. Some parallel 

f 

plate reactors can be modified to introduce Ar at low pressure 
and effectively as in a back sputtering mode. This latter mode 
could then be used as a clean up step prior to second metal 
deposition - although not in situ. Undoubtedly, dry processing 
is going to be extremely important in future double layer metal 
development work. 
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III. TWO-DIMENSIONAL MOSFET SIMULATION PROGRAM 
A. Long Term Objective 

* 

: )' 

The basic long term objective is to develop a two-dimensional 
computer program which derives the channel current in steady-state 
from given values of potential;: source, drain, gate, and substrate. 
Also given are physical parameters such as gate length, dielectric 
thickness, doping levels, etc, The substrate doping level should 
be variable in two dimensions and the substrate potential is 
independent of the source potential in order to Include body-bias 
effects, The program is intended to model MOSFETs with short Channels 
with particular emphasis on subthreshold and punch-through conduction. 

Secondary objective is to utilize techniques in the computer 
program to minimize computer run time while still maintaining reasona- 
ble simulation accuracy. 

B, Scope of Work 

Review the literature dealing with electrical models for the 
MOS transistors which takes into account two dimensional effects and 
investigate their applicability to short channel structures. Identify 
and demonstrate computer methods suitable for analysis and investigate 
preliminary problems which need to be solved. 
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C. Introduction 


At the beginning of this project no programs Were available 
which satisfied the specified objectives. There has been consid- 
erable work in two-dimensional modelling of JFETs including the 
works of Kennedy and O'Brien [3.1] and Martin Reiser [3.2] . The 
MESFET has also received attention in two-dimensional simulation, 
e.g. , Barnes [3.3], Various modes of MOSFET operation have been 
simulated or modeled using two-dimensional techniques, e.g. 
Schroeder and Muller [3.4] (saturation), Barron [3.5] (low level 
currents), Armstrong e al [3.6] (pinch-off), Kennedy [3.7] (effects 
of ionising radiation), Kennedy and Marley [3.8] (saturation). 

Mock [3.9] , Motta [3.10] , El-Mansy and Boothroyd [3.11] . Although 
these works all contribute to the knowledge required in reaching 
the objective of this task, only the effort by Mock is sufficiently 
general to meet the objectives' — and it is unavailable for public 
distribution. 

The work by Barnes [3.3] described the application of the 
finite-element method to the analysis of the MESFET. 

Further, the finite-element method has been used for some time 
in solving problems in mechanics and elasticity; however, it has 
only recently been applied to semiconduction problems. This method 
has the power to treat some problems, such as eigen-value problems, 
for which the finite-difference method is awkward if at all 
applicable. It can also be applied to the solution of field I 
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distributions governed by partial differential equations, and one 
of the most attractive features as compared to the finite-element 
method is purported to be the ease of treating non-rectangular 
geometries and irregular boundaries. For example, the geometry 
of the VMOS structure could be accomodated. It was decided to 
further Investigate this technique, 

The details of this evaluation are included in the appendix. 
In summary, the evaluation resulted in some skepticism that the 
finite element method would be sufficiently effective for semi- 
conductor problems to justify the effort, A more recent paper has 
reinforced this attitude. This paper indicates that the proper 

formulations of the semiconductor problem for the finite element 

* 

approach remain to be demonstrated, and, in agreement with out 
observations, point out that the application of Galer kin's method 
is subject to skepticism. Therefore, it was concluded that the 
further work should be based upon the finite-difference method 
which we have used before although the finite-element method is 
intriguing and may be further developed in the future. 

Key in the development of this task is the proper selection 
of assumptions. Certain assumptions are required to focus the 
effort and facilitate convergence expeditiously and, therefore, 
minimize computer time. On the other hand, too many assumptions 
lead to loss in utility. 

The basic assumptions follow: 

1) Recombination/generation is negligible. 

2) Hole conduction is negligible. 
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3) The conventional expression for Poisson's equation is 

}i . .. ' 

valid [see next section] . 

4) The mobile majority carriers are Included in an 
infinitesimally thick layer of charge at the Si-SiOj interface. 

5) The current flow in the channel is then described by a 
one-dimensional equation, 

. ... t .... 

3S __ 

J * q D -r- + q u SE 

n M 3x H n x 

where 

S » density of majority carriers in the channel, 

D - majority carrier diffusion constant, 

U “ channel carrier mobility, and 
E x ■ tangential electric field component. 

The first three sections are common in the stated literature 
and restrict the model from covering phenomena associated with 
breakdown and bulk generated leadage currents. The last assumption 
is key and separates this work seriously from that of Mock although 
Kennedy [3.7] used a similar approach for the sake of simplification. 

Key also are the normal assumptions which are not made and 
which are commonly utilized: 

1) the gradual channel approximation, 

2) constant substrate doping, and 

3) zero substrate bias. 

Although a constant effective channel mobility is utilized 
during the initial development of the program, it is not intended 
to remain a basic assumption. 
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Assuming step junctions at the source (or drain), the potential 
in the depletion region and oxide is defined by Poisson's equation 
with appropriate boundary conditions (Figure 9 ): 

1) fixed potential along the source (or drain) boundary 
(BCD and 1JK) , 

2) fixed potential deep within the bulk substrate (AL), 

3) zero horizontal component of electric field under the 
source and drain and far from the channel (AB and KL), 

A) zero horizontal component of electric field in the oxide 
far from the channel (EF and KH) and 

5) constant potential on the gate electrode (FG). 

For convenience the substrate in the basic program is defined 

* 

as having zero potential. The source, drain, and gate potentials 
are their respective applied potentials plus their work functions 
With respect to the substrate* i.e., 


f 

• V 

+ V„ . . 

source 

S 

Bi * 

¥ 

drain 

“ V D 

+ V Bi • 


» V„ 

Hh A i 

gate 

G 

mus 


For the NMOS example the potential 4' is always non-negative. 

The gate dielectric is assumed to be ideal with infinite 

resistivity. The charge density in the dielectric is assumed to 

be zero except for a sheet charge representing the silicon-silicon 

dioxide interface charge Q which is arbitrarily assumed to be 

88 

© 

50 A within the oxide at one grid spacing above the channel grid, 
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AB - YDIST - YDIF 

CD - JI *» YDIF 

EF - HG « TOX 

BC - JK - XDIF 

AL » XDIST 
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IKS = IK-5, IL5 = IL +5 R | 

JL<J> - JL-1, JK - JK1-1 L * £ 

JHAX1 - JMAX +1 
JMAXtj) - JMAX-1 

Figure 9. MOSFET Gross- Sect ion 


SUBSTRATE 

CO-ORDINATES 
A - (1,1) 

B * (1 , JK1) 

C - (IU, JKI) 

D * (IK, JMAX) 

F - (IK5, JL) 

G « (ILS.JL) 

I = (IL, JMAX) 

K W (IMAX*JK1) 

L A (IMAX.l) 
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Key to the simulation of short channel devices is solving the 
two-dimensional Poisson's equation without making unnecessary 
approximations. The finite-difference formulation for Poisson's 
equation is given in the next section. The majority carrier charge 
is neglected in this formulation except for the sheet charge 
representation at the oxide interface. The intent is to ignore the 
charge where it is negligible compared with the ionized doping 
charge — which is a good assumption everywhere except in the channel. 

One can then solve Poisson's equation for the two-dimensional 
potential distribution more easily and faster than if the coupled 
majority carrier current equation is included and the majority 
carrier density variable is included for every mesh node. 

After finding the potential distribution one can, after the 
fact, explicitly solve for the majority carrier densities and 
current in order to predict punch-thru and sub threshold conduction 
behavior. This approach is aimed at minimizing computation time 
while retaining the essential relationships to maintain sufficient 
accuracy for modeling purposes. One should expect a compromize in 
accuracy, under certain conditions, when compared with the more 
general approach of solving the two-dimensional conduction equation. 

D. The Poisson Equation, Difference Form 
The Poisson equation in different form is 

V 2 f ■ - P/e (10) 
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where 




V » electrostatic potential* 

" charge deneity, and 
e » permittivity o£ the material. 

For acceptor doped material the charge density is 

P ■ q(p-nN Q ) , depletion, (11) 

where 

q « charge density, 
p ■ hole density, 
n " electron density, and 
N a - ionized acceptor density. 

the hole density is defined by 


where 


P ■ N u exp(-u) 

u “ qf/kt , the normalized 
potential, and 


C 12) 


u * 0 in the charge neutral bulk. 

Consequently, in the bulk substrate and in the depletion layer where 
the electron density can be neglected, 

P “ “ q N q [1-exp (-u)] . (13 ) 

In the channel region the assumption is made that the electrons 
exists only in a sheet form of infinites tirnal thickness at the 
silicon dioxide-silicon interface. Since the hole density is 

insignificant at inversion, the charge density in the channel is 

Ps ■ -q(N a s+n s ) , channel, (14) 
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where 


N qS is the effective sheet charge due to the ionized 
acceptors within the grid spacing defined by the interface 
grid, and n g is the sheet charge of electrons in the 
channel. 

Near the oxide-silicon interface there are positive interface 
charges of density N ss » consequently, 


P m q Ngg * interface, 


(15) 


with the location being assumed to be 50 A inside of the oxide, 
Otherwise in the oxide the charge is assumed neutral, 


P » 0 , oxide, 


To solve the Poisson equation numerically, the 
can be expanded in difference form as 


(16) 
hand side 


V 2 , - V T V 2 U 


- V T "i.J-l + * "i-l.J + ° u i,j 


+ » u i+l,i t-8 u t>j+x ] 


(17) 


where V T » 


,“1 


A - 2 t»j-i(*> j _l + tyl . 


"I 


B - 2[W (W + W )] , 

JL jJU: ■ X X X 
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B » hj)!” 1 i and 

C * -(A+B + D + E) 

with the grid ayatcm defined in Figure 10. 

Using the charge equation (13) the complete PoisBon equation 
in finite-difference form is written for node (l,j) as 


.m 




qN. 


♦W'fV -T^ (18) 

The above nonlinear equation can be linearized by defining a function 
F such that 


F<U i,J-l» U i-l,r U i,j’ U i+l,j’ U i,j+1 ) ” %i 


w A 


(i»j) U i,j“l + B (i.j) U i _ 1J + C(i,j) U 1>4 


+ D (i,i) + E <^> «l,J+i 

- (qN a /cV T ) [1-exp (-U 1>;j )] - 0 


(19 ) 


Defining a set of reference value of U's and F such that 


Uj . . , U? . U* , etc. and Fj . , one can expand F around the 
x$ 3 ■"* J» x*^j* > j X 1 3 Xf 3 


reference potentials; 


r^+1 ■ p.r . 3F 

hd Fj " 


AU. 


SF 




AU 


i“l» j 


( 20 ) 


3F 


3U 


l.j 


AU. , + 


3F 


1,1 


AU 


3F 


3 U i, j+ l 


AU 


M+l 
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Figure 10. Grid system definition 
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where 


AU 


_ „r+l ,.r 

* Uj J , - u. 


ij-l iJ-1 w iJ-l 


etc. 


The above expressions defines a system of equations in A U as 


F 1,J " A« r J)4 0 1(j ^ + 5(i,J) + CU.J) AU^ (21) 

+ P«,J)4» 1W(J +E(1,J) 4U 1>J+1 + Pj (j 


where 


A " A, B “ B, D « D, E « E, and 


G - C - (qNa/kT) exp (-u' ,) 


i*j‘ 


with A through E being defined inequation ( 17 )„ The function- F is 

„r+l 


theoretically zero; therefore, by setting F . . ** 0 , a system of 

J 

linear equations in A U is obtained. For the node i,j 


Alu i.j-i + B4U i-i,j + clu i.j 


(22 ) 


+ D 4 U i+l > j +E 4 U l,J +l -- F 1> J 


XT 

where F. . is defined in equation ( 19 ) . 

•T»J 

By applying the method of lines, the linearized equations are 
resolved into a tridiagonal matrix form as 


A A U. , . + CAU. , + E AU . ... 
i,j-l i»J i.J+1 




where AU are taken as the most recent updated values of A U 
each iteration, but still keeping the U's fixed until convergence 
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is achieved. The above equation ( 23 ) retains the old values 
of AU along the x-axis and calculates the new values of AU 
along the y-axis . Obviously, equation (23 ) can be written for 
the inverse case. By alternating between the vertical and 
horizontal forms of the tridiagonal matrix, the convergence can 
be expedited. 

Special consideration must be given to PolsBon's equation at 
the silicon-silicon dioxide interface. Figure (HA) illustrates' the 
channel sheet charge S in charge per unit area, the acceptor 
charge N in charge per unit volume, and the discontinuity in 

£i 

permittivity. Applying Gauss’s law in one dimension with two 
surfaces at the midway points between C & C + 1 and C - 1 & C 
leads to the expression 


e E- In 
ox 1/2 


~ C Si^-l/2 


positive charge enclosed 


(24) 


N h . 

“ - q(S + } 


Multiplying through by ~ v ~ (h T h — 


leads to 


Si 


h +h . 
c c-1 


fox ( U c+1- U c) , U c-!- U c 

e Si h c 


c-1 


la. 




N a h c-1> 


r V 

On j T ft 


'' S i T 


h . h 
c + c-1 


(25) 


which is the one dimension form of Poisson's Law 


dy 


eU- 


peff 


e Si J £ Si 


(26) 
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with 


2q(S + N h /2) 


Therefore, the finite difference form of the equation is obtained 
by modifying the following terms in equation (18): 


E * 2 


/ \ 
e ox 


"si * 


iVVi + V r * 


( 27 ) 


C » -(A + B + D + E) , and 


P eff _ 2q (S + M a h c-l /2) 
e Si “ Vsi " <V h c-l> 


[.Note: y ■ J for the channel node in the program.] 

JTtdX 


Similarly the linearized version is modified by 


h m V \ m c i» and 


(28) 


F i,c - Vi,c-1 + B i U i-l,c + C i U i,c 


+ BjU. . , ^ + E 4 u, „ x ,+ P -//e. 


i u i+l,c "i w i,c+l' eff '“Si 


where the subscript i has been Inserted to denote the coefficients 
for mesh node (i,c) . 

In the region of Si - SiQ-j interface where the Q gt , charges 
are assumed to exists, one mesh node is assumed inside the oxide, 
the coefficients are* 


A « A, B = B, C 13 C, D * D, E .* E, and 


39 






r - Au. + BU 


.+ CIL DU 


+ EU. 


i,c T DU i-l,c+l UU 1 > C+1’ i+l,c+l l,c+2 


2Qgc _l 

+ rf“ [h c + h c+i ] 
e ox V T c C+1 


( 29 ) 


Similarly, in the oxide where p “ 0 , the coefficients are 


A - A, B * B, C “ C, D “ D, E - E, and 


(30) 


F ■ + + «i,j + Du i + lj + Eu i. 


E. The Channel Conduction Equation, 
Difference Form 


As noted previously, a significant assumption utilized to* 
facilitate convergence, is that the channel charge is modeled as 
a sheet charge at, the silicon-silicon dioxide interface. The 
channel current is assumed to obey the single dimensional form of 
the conventional three dimensional equation, l.e. , 


n 


q D — - + q u S E 
4 n 3x 4 1 n x 


(31) 


where 


S « 


D 


n 


n 


E 


density of electrons in channel (//per unit area) , 
electron diffusion constant, 

electron mobility, normally field dependent, and 
tangential electric field component. 


The above equation can be written in normalized form utilizing 


\ * V T lx and ^n 0 ^T^n (Einstein's relationship): 



«*c» 


- 9 S __ DUv 

J m Q V m U (’T — - 3 *r— ) 

n 4 T n v 3x 3x' 


fg " (S-US)M 


where f 


V" V" 1 - s '- H ■ “" d «'* i 


At any point x in the channel the general solution of S is 


S(x) ■ exp 


L u’dx | / ( -j— exp / u'dx) dx + c| 

x i-l ( Vl n x i-l J 


which reduces to the specific solution for S(x) , given 


S at x » x, 


S( X )0-“ . ^ + A-. J x 


-u(5) 


Rewriting the above leads to in terms of : 


S 1 exp (-tfj) - S i-#1 exp <-U ± _ 1 ) 


M l i 

exp|~U(€) I d £ _ n 


x i-l 


Similarly, integrating between x i and X^ +1 leads to 


S i+l exp ~ j U i+1 ) " S 1 GXp ~ [ V - , J n 

r x 44.i m 


exp-U(S) d g 

y(«) 


Assuming current continuity, i.e. , « constant, leads to 




(38) 



***‘.±1 ScA -> -. v. 


i+1 ) 

-* wj , 


- <exp 


+ < exp - U A 


5 1+1 


EiElHiS A £ 
y(5) 


i-1 




- <exp“U 


1+1 


‘i «^Q|d« U*.-o. 


1-1 


The above equation is of the form 


a i S i-l + b l S l + C l S i+l " 0 


and define a tridiagonal linear system of equations. 

Equation 36 should reduce to a familiar form as x^-x^j" Ax + 0> . 

Let S. - S , S.“S + AS, U. U , U « U +AU . Then 36 becomes 
i-1 o i o i-1 o i o 


-u -u 

(S o + AS)e °(l-AU)-S o e 0 


-1 -U Q 
M e 0 


Ax 


qv« 


C 39 ) 


which reduces to 


J - qV JJ i~ S - S ® 
n M T \Ax o Ax } 


(40 ) 


which displays the familiar diffusion and drift terms in the limit. 

■ i ». 

The integral form in (36 ) has the potential of being more accurate 
with larger grid spacing. 

Another observation one may note is that the coefficients 
a l * b i * c i ma ? ^e sca l e ^, or normalized, to values less dependent 
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on the exponential of potential by multiplying them by exp (2U^) 
which gives 


rx 


a i " " exp(U 1 -U i-1 ) I ^(OexpU^-lKOJdC , 


v l+l 


M" 1 <5)exp [U,-U(C)J dK , and 


c l-l 


V“ exp(U i“ U i+1> 


1 il*~ 1 (C)exp [U^~ II (5)1 d § . (41) 

‘i-1 


And, of course, if the mobility Is assumed constant, then it can be 

removed by multiplying through by p . 

Noting that the Integrand varies exponentially with potential 

whereas the potential Is expected to vary with x no worse than 

the third power, one can approximate the "integral by assuming 


U - U ■ a(x - x ) 
o • o 


dU 


(42) 


and evaluate the slope a » ^ at the lowest value of potential 
over the range of integration. In this manner the neglected higher 
order terms contribute very little error since the integrand falls 
off exponentially with potential. Changing variables for the 
integration of b. , as an example, leads to 


where 


1+1 -1 


U 


1-1 


dU 

dx 


P exp [U . - U] dU 
& 1 


Further, 
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(44) 


b A " [apT 1 {exp (Uj- U i-tl ) - ex P^ U i” U i+i>^ 

with a and p evaluated at the lowest value of U which for a > 0 
is . Similarly, 

a i " l “M 1 1 ex P ex P( U £" %.j) 

-3 

c i - tap 1^11 - expCu^ U i-X ) 1 exp(U i - U i+1 > 

-I 

b i " l a PJi_i( exp ( U i” u i_i> ~ ex P( U i" u i_i> 1 ( 45 ) 


where the subscript on [ap] denotes the node for evaluating a 
and p, assuming a > 0- . 

The boundary conditions for this system of equations 
association with electron conduction in the channel must be defined 
in terms of the source and drain corner meBh nodes for the MOSFET 
operation in the active region. Figure 11 Illustrates the applica- 
tion of Gauss's Law at the source corner of the channel. Mesh 
nodes along the boundary of the source are common pqential, therefore, 
the tangential electric field terms along the boundary of the source 
are zero. Since the electric field along the channel F,g for 
operations of interest - at least at this source corner, the surface 
charge is modelled in a single dimensional mannec, i.e. , 


e E M ® - q S 
ox N v 


(46) 


or, in terms of potential. 


Vox ( °r.c-H~ U r.c> 


* q S 


(47) 
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An alternate scheme Involves the solution of the one 


dimensional Poisson equation: 

where V Q « the applied bias to the gate 

■ the gate work function relative to the substrate, 

V c ■ the applied channel potential, 

« the work function of the drain/source, 

t „ ■ oxide thickness, 
ox 


“ %<,- I s + q, < 48 > 


1~ I <V ♦»> - <V V 


e “ oxide permitivity, 

OA 



silicon dioxide-silicon interface charge, 


q ■ electronic charge, and 
Qg ■ the bulk charge density (per unit area) 
■ 0 for drain and source. 


The Qg expression can be defined in terms of the vertical electric 
field in the channel region: 

Q = eV m ( u c-l ~ u c) - q Na jigzl 
■B T h , 2 

c-1 

where the subscript c refers to the channel node. 

F. Channel Coupled Equations 

Although the channel conduction equation .(38 ) is linear in S, 
it is nonlinear in potential. The simultaneous solution of the 
conduction equation and the Poisson equation along the channel 
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requires linearizing (38). 


Ut F' t - S w + \ S t + =i S 1+1 wh.« a v b ± , c x 
are defined in (45). Since the terms [op] vary more slowly 


than the exponential terms , the term is assumed constant for any 


particular iteration, using the last known values of u*s for its 


evaluation. Then, the linearized conduction equation becomes 


0F r 

A ci “ 3^ ■ ■ ‘i s i-r tou, i-i s i o * p<u r u i-i ) 




B ct'“t > 


°ct ' S& ■ - Vi-i' oxp <VW- l i + b i s i 


f Vi+l 11 ~ c *^i-l» l> 2 expiO-U^)), 


D ci " b i * 


1+1 


+ bxptUj.-U^) (expCU^.p - OKPCVWr 1 


" C i S l+l * 


G ci “ c i ’ and 


C- F cl + A cl 4B l-l + B =1 4S 1-X + C cl 4B i + B ci 4S i «» 


+ E ci AU i+l + G ci AU l+l " ° * 






The Poisson equation has been previously derived, equations 
(27) and (28), with the coefficients being defined below for t 
new system of equations: 


where 


B i 4 u l-l,c + V u ic + "i iS l + Vw,e 


•V u i,c-r ¥ u i,=+i - r lc ■ *L 


(50 ) 


-1 


"i ' **" ^X.O-X^ * 


A 1 " A i from ( 17 ) » 

from (17 ) , 

D -D. from (17) , 

1 

E. » E. from (17 ) 

1 X ' * ' 

Ci - C, from (17 ) , 

F. r « as defined In (28 ) . 

X j y C 

Simultaneously solving for the channel potential and the 
change charge requires the simultaneous solution of the above two 
equations (49 ) and (50 ) which results in the following system 
of equation in matrix form: 


C H D 0 0 0 . . . 
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The matrix Is a 6 diagonal matrix, The number of variables Is 
twice the number of grid points within the channel— not counting 
the boundary nodes at the sourco or drain. The solution can be 
conveniently solved by using a Guess Ion elimination technique, 

G*. Overall Flow Chart 

In order to expedite convergence, Initially, the gradual 
channel approximation Is assumed and the potential Is calculated 
along the channel, This analytical approach divides the device 
cross-section into two regions, the channel, silicon subBtrate 
region and the oxide region. 

In each of the two regions an initial approximation of the 
potential distribution is obtained bv solving the one-dimensional 
form of Poisson’s equation utilizing the given boundary potentials 
including the channel potential derived in the previous step. The 
one-dimensional equation is equivalent to setting dU/dx “ 0 . 

The potential distribution is found along one column (vertical) 
simultaneously. By sweeping in the x direction, i. e., repeat 
the solution on another column, the potential distribution for 
the entire region(s) is obtained. 

Next the solution to the two dimensional p, E. is obtained 
in each region. In the oxide the potential along a column is 
obtained with the region being swept in the x-dlrectlon. In the 
silicon convergence As expedited by using the alternating direction 
implicit 4^1) scheme. First, the potential in the row below the 
channel (JMAXO) As obtained. Next the potential distributions for 
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each row between the source and drain are obtained by sweeping in 
the negative y-directlon. When the solution to row JK1 is obtained, 
the solution to the column I n 1 is obtained The colump solutions 
are obtained by sweeping in the positive x-directlon. This ADI 
scheme is repeated until satisfactory convergence is obtained. 

With the two-dimensional potential distribution obtained for 
both regions based on the gradual channel assumption, a new 
estimate of the channel potential and the channel charge is obtained 
by solving simultaneously the P. E. and current equation along the 
channel row. The 2-D P. E. is then solved along the rows between 
the source and drain by sweeping y in the negative direction. Next 
the columns solutions are obtained sweeping the x-dlrection for 
the entire regions — thus crossing the channel in the appropriate 
locations (I » IK + 1 to IL-1) . This ADI scheme starting at the 
simultaneous solution of the P. E. and the current equation along 
the channel row is repeated until satisfactory convergence is 
obtained. 

The channel current is obtained from an explicit solution to 
the channel current equation after the channel charge and potential 
are found. 

H • Progress Report 

' ■ .... • , 

The basic program has been written which provides user 
flexibility for defining the physical dimensions, the number of 
grid points, p or n type device, metal or silicon gate, type 
orientation of substrate, etc. Provisions have been incorporated 
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for sweeping multiple bias points in order to plot, for example, 
drain or transfer characteristics. In addition, provisions are 
provided to plot the two dimension potential distribution as 
well as, print out this data in table form. 

Using an example bias the program has successfully completed 
the steps through the solution of the two dimensional P. E. using 
the fixed channel potential as a boundary condition. However, at 
the next step, the channel charge is found to diverge. Time has 
not allowed the resolution of the problem; however , an error has 
been found in the current equation program which may be the cause. 
Of deeper concern is the possibility that the channel charge is 
not directly tied to the channel potential through the P. E. , but 
only indirectly. This may: xequirs; ■..jfch.fe formulation of an expression 
for channel charge qS as a function of potential and electric 
field which Is then inserted into the PoiBson equation. Further 
investigation of this problem is required before the best solution 
can be determined. 

!. Availability of 2-D MOSFET Simulators 

The first 2-D simulator program which could provide I-V 
characteristics was developed by Mock [3.9] and Kennedy [ 3.8] at 
IBM in the 1970 to 1972 timeframe. The program, based on finite 
difference techniques, is proprietary to IBM. However, a modified 
Version of this program has been developed by Sutherland ['3.14] 

(and Kennedy) at the University of Florida under contract to ARPA. 
This program is restricted to uniformily doped substrates and zero 


substrate bias. 

Kennedy at the Applied Electronics Research Corporation has 
modified the Sutherland program to Include provisions for 

1) 1-D substrate impurity profile, 

2) mobility dependence on gate voltage, and 

3) substrate bias 

and is now licensing the program (NEMOS) for $10,000 for one year 
and $7,500 for each succeeding year. 

In addition, Hitachi has developed a finite element version 
with provisions, at least, for substrates with 1-D Impurity profiles. 
Dr. Asai of Hitachi has Indicated the possibility of licensing the 
program to Universities, government and industry in the United 
States. The Hitachi program called CADDETwas developed in 1977-1978. 

At least one other 2-D simulation program is under development. 
Jim Melndl at Stanford is under contract with ARPA to continue the 
work by Sutherland. 

At this writing the author has been unable to determine when 
the program under development at Stanford will be available for 
distribution. 

J. Recommendations 

In view of the recent "availability" of alternate sources of 
the 2-D simulator, it seems prudent to evaluate more seriously the 
availability and cost for achieving a satisfactory and working 
program from the four possible sources. The compatibility with 
existing computers and the difficulty in modifying them to meet 
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the stated objectives must be considered. The authors believe 

ii . :s 

this should be accomplished before seriously continuing with the 
present program. After all, the 2-D MOSFET simulator Is just the 
means With which to evaluate various short channel MOSFET device 
structures and not the end product. 
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»3Pc2rLoif iitixliu 

CALL PLOTulpcJoIliitiO* ■ . 

S ! h-<s:i3s»oi . 

DO 17 J-ti ICtNYSTP , 

so vo >oo 

■0. I A I TSPf l , «0B » VAL ,n. »" 1 1 

*°° CAbL N N{|HBcKi-oU4,YSPc{iO*|Orl0.oO*i»|1 

IT CoJ»T|NUC ,E * , " 0 * SS * VSfC, *' 0, ' 0 * , ' ¥ * , ’ ,,0 "* n 
.. WALlivALI*tO. 

I A jKa JK* I 0 

XVaXHAXV/FLOAT(lG*l0l 

|6R103>IGRIn*l 

!8>« l |:||lS5l85„.s„ 

XSPCI-XSPC-0.2 
. VAL "XV • I 1-1 I < 

I* CALL »!UMHERtXSPCI ,”0, >0,0*1 ,VaL. ft. 0,2| 

••• PUT LABELS 

CALL SYMBOL <-0«A, | »0.n* 131 3. IV, 90. ,30) 

CALL SYMBOL! 2.00,»0«Sn,0. I 313, 1 X .0. ,30) 

CALL SVMBOL<2*CoS,S.Hfif ,• TEMP ■ 1.0, ,71 

CALL NUMBEfl 3,0,S.0,Hr.Y ,TCMP ,0, , II 

CALL SYMBOL I 2, 0 , 5 , 3 , hc»Y » * MOBMxL 1 7 1 N6 TIMr • 

• 0.0.321 , .. 

CALL NUMBCH(3.A,s.3,h«T,TM*X.0..|I 

C*LL 5YHB0L(2.5.S.| , N«T,.* M ji t NT a »,0..lO> 

CALL SYMBOL(3,O.S a 1,HloT,|L.O.,l7t 

IF (N.EQ.2) CALL S THRO!. « 4 .1 , S.i »MGT , * TmI C vMCSS • 

CALL NI|MBebIS.A.S.B|M6T,VI)|0T»0«« 1 1 
••• ORAN CURVES 

vlogbVMax /float < j c R i d 
xvi»xmax/xm*xv 
VMVR al ./YM|N V 

00 22 I ■ I , N P L T 

CALL PLOTIO.O.O.O.3) 

IFfl.EQ.il NST P»NS JP I 
IFH.EQ.2l nSTP»NSTP 2 
IFll.EQ.3l HSTPaNRTPO 
IFH.CQ.H) NSTPaNSTPM 

1 Ft I . F Q • S I NSTPbnsYPS 
00 21 J»I,NSTP 
IFIJ.EQ.il L"3 
IFIJ.NE.il L»2 u , 

IF I IAMBNT.NE.3I T|Hf H*JI»T1NC! I , Jl/T |MC IMPLY, NSYPl 
XMOVEa TIME! I .JltxVI 
IFIN.EQ.3I GO TO >01 
IF ( N . F Q « 2 | GO TO 702 , 

YMOVEaYLOG'ALOGIn, RS( , , JlavMVR I 

201 TMOVEa?LO<i»*LOfi I 0 I Q ( I J I • Y M V R| 

, GO TO 21 

202 VMOVCaVLOG«AkOG|a|Xj< | .JiaVMVRI 
21 CALL PLQTIXmOVE.VMOVE.L.OI 

IFIN.NE.2KG0 TO >2 
CALL PLOT 1 0,0 ,0*0,3 I 
00 23 Jal.NSTP 
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triJ'CQ'll L»1 

If! J.NC.I) C»2 


IMOVC>T|Hej I 

ri’SiriWBMlU 

it g^^aMBKWjtar*"' 


1 1 1 « ji • i »r»oi 


. cut PLOT <0* iO* »t0?l 

M PONT |NU{ 

«TOP 


PNO 


( 
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S i 1 


l|C»ON*SnS?i I I.PAHtM 


i! 

13 

1*1 

Ift 


V 

iCt 


SUBROUTINE P ARAM I f I i T I 

ILLICIT DQUBLC Mtc'slON *A*HtO*Z) 

Efi>i*n» 7 .ic»iOMART<rt) 

fG»CG/t«» A2C-6* (T I I 

awsHu 

ir .ANft'.R.ST.I.OoH «0 TO I 

RETURN 

FNO _ . 


s 

8 

8 

0 



1|CR0N«SflS2l I » iC 


DOUBLr PRECISION rUNctlON filfN* J.i) 
IMPLICIT DOUBLE PRCC I S I ON U-H <0-7 | 


IRRlN'S conductivity PORMULAS, 


i r i jj.eoti i go To, a 

iriCN.GT.o.m a» t ,0 
IPICN.GT.n.nl O-T.IO-tT^ 

IP ICN<GT • T • RD MM A.n<4S 
IP jcN<GT<l<SO ♦ia*?*3.3d-i» 
IT (CN.GT. J.HO*»A> A»0< *32 
iMCN.GT.2<‘lO<|ll»»l ,*70-|A 
IP (CN.6T<I<S0<IT) AbQ.TAA 
|F(CW.6T»l.S0*IT|»"**,r>“l7 
GO TO & , 

IFICN<GT<0<OI A;| ,0 
iriCN<GT<0<0)R>2.n-li , 
IPICNtGT • 3*50 |S I AaOt 0)7 
IFICN<6T<3<S0I S |0»A»f T0“| *1 
IPICN<GT<I<0DI7I «bO*gH 3 
|P(CN<GT< I < 00 1 7 IOa4<yjD-T 
IPICN.GT < Y « 60 I • ) *“0<f *♦ 

|P I CN<GT<T<SO I •|«»2<0- 1 A 
IPICN<GT<A< 0 1 T ) »«0,?H*t 
IPICN<GT<A< 0 1 T )Na I <NjO- | 2 
IP ( CN<GT<2<3S020 1 *»0<NiA 
IPICN<6T<2<3S020I •■I.ONO-A 
« m R*ICN**A> 

RETURN 

END 
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4|C*ON*SnS21 » I .0*0 aTA 


I 

r 

* 

i 

i| 

ii 

1 ! 

K 

it 

ii 

it 


SUBROUTINE (iXOATa|AHrnT,ONINT.T,h,C.H»KB| . / .... 

SUftftOUTlNC ron OXinAtlON PA(U»C?t*$ AWO SC«"C§a?| 0« co*frrift 


IMPLICIT OOUOLC MCCtSlOHU-H.O-n 
intcscn ahbnt.okinT , 

DOUBLE PRECISION *B #M ,M | .M3 


tr4AHHNT.tO.il SO 

MN.ROfTTo-lOaPCXR 


I; 7 * 35 . 7 *|/TI 

Sfi!:i«:!i:ii i$K !««»• 

60 TO M 

f;Kinbs:j?‘ 8 si:iBiwsKtfi..- 

!f!S!!!l?:f 8 : 5 i I:!:S« 5 {{:|: 8 B«:I 5 IH:JKW. 

CONTINUE . ...... 

MI»33.3»DCXPl-0.S»/ KR*TII 
N3*20.0*0e>P(-0.67/(K6«Tll 
iriORINT.C0.il NrN* ..... 

iriONlNT.CO.ZI N«|NI«*3!/2.0 
jr (omnt.e0.3i h.h3 

CNO 


If 


IN 
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R|CR0N*Sf)S2l 1 1 


2 

I 

R 

I 


5 

ift 

1! 

if! 

22 ! 

« 

a 

n 

n 

U 

31 

ll 

3 ? 


10 

IX 

« 

#3 

11 

1 ft 

14 

17 

1 ft 

If 

*o 

si 

13 


•FRONT 

SUBROUTINE FRONT! |MAX| .JHAXI ,l» ,X , TEMF ,DF I ,DELT , JsTEP , 
• KTTFE.KOT|.XO|StI 

gun. FOR CALCULATING. CONTOUR FRONT MOVEMENT DAT* AND 

Hading and writ in* tmf sane on file 

INFLIC IT DOUDLC FrCC 1 G I ON I *»H , 0«f | 

01 MENTION AFRONTIaJ* TfRONTlAI 
DATA I DOL , I aLN* , I AR / 1 H« , IN ,lHc/ 

IFlXTYPE.NC.il 60 TO rOO 

read file 

00 AO KKsIiJSTEF „ 

READ IY.JrOIMArK 
if imarx.ne. idol • go to to 

!ES8I5:558I WK"' 1 *"'" 

READ IT. TSOI lOUN.Lal.AI 
READ “ “ ‘ * 


SO 

AO 


REAO I f « T26?*?HAX O.UMAX I *K .I.M 


INUE 


RETURN 
•00 CONTINUE 

00 ZAO LLal.A 

10.0 


2 IP 

240 


270 


T 20 

T30 


CONVACal0.O*«FLoATl ?0-LL*l I 
XFaCONDEPl IMMl ,JNAil .O.-UMAXl .O.O.CONvALl 
TFaCONOEPI IMAXI ,JMAX| . -2 .O.O.O.CONVAL I 

IF (XF. LI). 0.01 rO To . 250 . . . , , 

XFRONTILL|a|XF-FLOAT«XI I • I KO I ST/FLOAT I |MAXl"l I | 
CONTINUE 

IF (TF.EQ.O.OI 65 Tn 2A0 

TFRONT<LLI«IFLOaT( JnAX| l-TFI/FI OATI JNAXI-I I 
CONTINUE 

STORE CONTOUR FRONT MOVEMENT DATA IF IFILr a | 

IF*UO??»«T. 1 1 MARX-IrLNX 
WRITE I V . 3?0 I MARX. 

IF IK0TI.6T.I I 60, TO f7o 

WR I TE I T . T30 1 TEMF ,DF I ,DELT 

MARKalBLNK 

WRITEIT.JTOI MARX 

WRITE IT.JJOI IXFRONTlLLI.LLai.Al 

WRITE <T» NSO I I V.F'R 0 N T ( L L I . L L ■ I . A I 

WR I Tr I T . T20 I IMAXI (UMAX | ,K.LM 

FORMAT I I HO • Ail 

FORMAT I I M .3Fl0.il 

FORMAT I IM « T I S I 

FORMAT I IN .AIEIR.T.2EI | 

RETURN 

END 


^SStSSfJ! 


166 


1 1 CttON*Sp52 1 1 ) »FRRLOT 

•••PLOT f«06H A M FOR TH t CONTOUR FR, NT HOVEhENT ••• 


f : i 

< , ;> 


I 

i 

R i 

ft 

Ai 

l 

10 

II; 

I! 

ift 

!! 

to 


n 

H 

U 

It 

i? 

32 

33 
3 H 

n 

37 

31 

3t 

» 

tl 

*3 

ts 

t4 

tr 

H* 

tt 

so 

u 

53 

s*< 

ss 

54 
ST 

St 

St 

s 

43 

44 

45 
44 
47 
44 
4t 

71 

7§ 
74 
7 S 

74 
7 7 
74 
7 1 


C 

i 

i 

c 

$ 


RE An I N FROM nATA OECK 
jSTf f • « OF TIME rTEPR TO AC READ IN 

♦TCHR - SITUATION TEMPERATURE*, 

♦Off - » AHOA*«t«i 
*OLT - TINE STEP • » 

OOUaLC SJCClfigN iFR0UT<*00»A| *TFRONTM0n.ft» 

double precision tchp.ofj.olt 

OIHCNSION TaU(«0(|I ! 

OI9ft N ?J9.N JKJtJilVPi.jAiftl 


IDOL, Ihi 1 n|I ^/{hSiIm *lH»/ 

IaFAHCONcCN, AHT R amO i AHN F RON $ AH? HftV| ,AHHCNY P,AHLdT 


OIHCNSION 

oata 

! ..o IMMImU 

RIO rORHATUH *R|ftl . , 

R90 FORHAT ( I H iAIClR.RttXI I 
C •••READ DATA FROR nCCK 
, READ TOO , JSTER 
200 rORHATISI 101 
TIHE-O* 

RE AD I i • 3 *H ARK 

lEiSK^SSi'JSW.SliMtl 

4 slisisisisi !?flssUK:tt{:tt:i:t! 

RE AD I A «R2Q I IHAKl , JH a«I |UI 
yihc»timc*oi t 
. T AU I RK I »T I He 

1 CONTINUE . .. . . 

OATA HGTtXMAX,YHAl/0.n*7S,7.0,ft.n/ 

C ••• INITIATE THE Rl,ftT 
00 ID N-1,2 
CALL PLOTS! I 0 « 0 , I f|«0 I 
, . CALL PLOT! I.S.I.R.-S) 

C ••• DRAW BORDER 

CALL PLOT (O.O.YHAX i 2 , • t I 
CALL PL0TIXMAXiYHA*i2,Til 
CALL RLOTIXMAXi0.ni2.*l » 

CALL PLOTIO,0,0.0,2»-| I 
C ••• SCALE GRIs) 

TTAUaTAUt JSTERl«lnO. 
lOlV-INT! TTAU) 

IS»I0|V/2S«2S*2S . 
TAUHaFLOATIISI/iriH* 

WlaXFRONT ( t til 
IFIN.E0.2I V Y I "YFrONT < I » A I 
00 2 | ■ 2 . J S T E P 
VY2-XFR0NT < | ,A) 

IEIN.E0.2) YY2«YE*0NT ( I t At 
. V YVpAMAX I I YY I * YY2 I 

2 YYI-YVV 
VAP«YVY»100. 
lYO-INTlYAP) „ 

|T«IY0/25*25*2S 
YOHAX-FLOAT! ITI/lrOt 

C ••• DRAW GRID 
OXpXHAX/20 • 

00 10 I ■ I • I 9 
IpFLOAT! I I *ox 
CALL plot I X ,0. , 3) , 

„ CALL RLOTIX.YHAX.Ttl » 

10 CONTINUE 
OY-VMAX/20. 

00 20 J>l • 19 
Y ■FLOAT! J»*OY 
CALL PLOT 10.0 ,Y ,31 
„ C ALL PLOT I XM AX • Y . ? 1 1 ) 

20 CONTINUE 

c ••• erase labi e arf* 

00 30 I ■3,1*1 
X X • ! I ■ 1 1*0 X 

CALL PLOT ( Xl , I 7 »*«Y ,3| 
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• 0 

SI 

u 

R 

ii 

8 

ii 

i! 

♦t 

|00 

IB 

1 07 

lot 

I Of 

ii! 

ii 
i 
i 
i 

!»» 


30 (ILL PLOTlXX.If.SiDY.fS) 

ia citt 

iA. Mil* ff U U k Ml 


C ••• PUT S THIOL 


rw| sinwv i . * 

CALL >VM|0L4 Q« ft ■««? »0* » *11 * I A .0.0 ,33 1 

iriN.re.n ixair , . 

CALL ST MOOLI I »A»f* *•*»«< 


•••AXCS NUMirOS 
V*- 0.2 


» .A* 13 1 3 . 1 Wt0*0» I ft 


.. WALafLOAtll- 
*0 CAL “ 




NUHOeoli tV*HAT«VAttOtO#l| 
VAL-FLOAT | J-l >*V ohAX/? 0 , 

.Vi i iillUaf ■ t W W ii 1 1 tl k 1 . t 


w-.iL N 


w.w, ANfC* 

NORMAL ! ICO TIHr • VAUNftttSl) 


40 Hit SVB*SCi::?*r*!:riilS:?*Ht 

CALL STM0OLII»»“«l»»|ll3#*N0OHA 

c2u 0 PLC!TUto.0.ljt3l 

cSll £ plotuau,xhovc,2;oi 
00 70 HL«1»JSTCP , 

TAUP»TAUIHL»*XHAX/TAUH m . „ „ 

IMOVr»XrAONTIML.LL> •YCaX/TOMA* 

IPIN.ee*>) VHOVe»rFgONT|ML,Ll»*m«/YOMAI . „ 

IPlTAU.fiT.O.O**NO.THOVe*eO.O*0«ANO»N«CO*»l CO TO 70 
iriN.ee. 2».XMoye # vHgwr- „ 

|plXM0VC«LT.0.r XN<>ye«0.C 
, CALL PLOTITAUPtXHnVE.PtOI 
70 CONTINUE. . 

CALL PLOTlO..O.»ftf » 

00 CONTINUE 
ATOP 
FND 
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.! .!!• I Ml MM!) 




l|C« 0 N«SoS 2 lM*SMr uc 


1 

2 

li 

l 

t 

t 

10 

tl 

u 


i 


RURROIJT I Nt SMr Ilf, if iHr# OCWT *0*1* Y » ***** f H« 1***1 v J I 'fjUHC •' 
• rJHAXl itt>lST#*ft*IAH»|Nf I 

SUr RRtTCS Seer RfrISTiNCE* JUNCTION OCRTn |Nr INTliRRTCO 
iRfURITY ON UNIT M 


l”WlllJonj L j ? "“u«°l!} 

m SMW5 Kltf l • . f H.« ..CCT . Otc » , Jgi *T 

RflTf ! ! ? • iOOr , R5| V J I • V JUNC *0 » TINE* *0 


1 13(100 


i 


100 foAnaTI IHO tRfl 
RETURN 
CNO 




I ). OUTPUT 


1 

2 

* 

S 

7 

t 

10 

ii 

IN 

il! 

1 } 

if 

1 ! 

ill 

27 

29 

If 

H 

)l 

2* 

27 

29 


y . 1 H A x | ,JMA*I ,K.i«. 4 j»nHt#Vot»T. 


c 

c 

c 


SUNMOUTIM OUTPUT I * I V • I H A* I • JM A 

• I 0 1 I T I NC » *0 * P T I Mf , 0 1 1 n € . I AM9NT I 

TMANSIENT OATA PH »nT Ot|T 
IMPLICIT DOUILC PtEC M!|0NU-M,0»n 

snsiMm 

pSISI 00 1 iMjtlMffPTlME.OTIMEiMilUl.l^Mi 

PM I NT (Oil INjNbJ, IMAK l.il 
PM I NT 1 02 • I X ( | ) . f *2 « f MAX I • 2 1 

Sibils 5stft!2i>' 

IP 1 1 *MPNT *Nf *2 1 Q»(JHA*I“J|*4 I • /KLOAT | JM*XI*| 1 1 
2 PM | NT 103. Oil CHI 1 1 I Jl I 1*2. IMAKl .2) 

. :»{*S f'Asigf- si: 1 ;?' , 

18* POMMifnM0H0X**»*0X |f>E TMICKNESS*IN*CM « * / 1 HO 1 1 3X » 1 1 1 1 PP 1 0# * » ) 

100 p ORH AT i?H6i| , 0X*l?yT|Mr5TEP ■ . | A , 2K . 7HT » ME • |P)B*X. ,,, 

•SX.»ClAPSE0 TIME »N SpC* I'.SX . 'PMCOCP ■ * »P 10. J . 2* i »0M I Vf IN ■ 

• Pi 0*27/ 

• )flX.»OXlOE POSITION'/ „ _ . .. .. «... 

•|0X.*X<*.|2.«> ■ . iCA . 2 . 2X . • X I • . 1 2 » • I • ME*»t// | 

1 0 I POMHAT ( // . IMO. JMI • .*X,I2|I0I „ . 

POMHAT ( /» lH0.2X.SHX • „ . 7 X . I 2 1 IPP I Q. 2 > I 
POMHAT I I H. 2X ,» Y-l • IPP*. I .2X , I I IIP? 10.21 » 

op TEMaTION > ••IS) 


* I 


18 ! 
1-0 M 
»os 
10* 


POMHAT I // . IMO. *N0« 

POMHAT U HOI ' . .. 

rOPMATI IHI 1 1 OX 1 1 S»M »TP*. 'TIME * . AA I 
MCTUMN 
ENO 


PMT SoS7.PiOT-CONTOUM..SU9IONi.TM|hAO..AHC..XV2..PLOT..COnOEP..MAIN 


*|C*0N*SnS2l II 

I 


• HOT-CONTOUR 

*UMOUT|Nt PLTCONjA.coNVAl.NC. In , jA I 
IROCONCCNTRATION (M. OT rU|RR0|RaH t 


Of MENS ION A( AS|AN| 
OIMCNS ION CS i 2 • 3 I f 
DATA llERO/l .OfloS 




»ni w » i 

S:?:8: 


i-O.I.O.O.O.*/ . _ . . 

j^9 • 0 1 • 1 1 0 • 0 • 0 • n • • • I • n « 

.•0»0*9iIiQiI*Q«0*0/ 


lUARE C*0SS|NSS 


,M« nytf 

•o*Of-{»o»o*Q«-(*o*o«o*!*o»i*o*o«n«-i«o»o«o/ _ . , 

Hit if i Ii{ itt: i: ' 

BATA It M A 31 1 VMAX/ 1 * m ^ § n/ 

C§n1uR- ALOG lO(CONyAL) 

SC ACX— XH AX* 1.6/1 1 V "1 1 

» C S V A R T ^ C 0 N T 0 5 1 * At ARCH 
NT*0 
It* IN-1 
. JL- JK- I 
10 00 §0 l-2.lt 

• UOCATC loUARC C*OSS|NiS 
11*1-2 

iV ■&• 26* I A tOft 1 0 1 1 1 1 1 41 f * AL06 i 0 1 A 1 1 * | • J » I * ALOC 1 0|A 1 1 , J.| II* 
*At06| Ol A 1 1 * I « J* 1 1 1 1 
Ml >-IO.O-*M 1 1 
MRI-all) _ , 

„ • locate TfiANHrs 

20 00 SO K- 1 » A 

NC-0 . 

.... 60 TO ( 2 1 .22.23 • 2S I *K 

21 R(2I-A( !♦! .ji 
M3I-AI I .JI 

_« 60 TO 30 

22 M2I-AI I • J I 
M3I-AU * J*H 

, 60 TO 30 

23 M2I-AI 1 1 J*l I 
Rl 3I-A< I ♦! » J*.| I 

u 60 TO 30 
2H 8 ( 2 I * A ( 1*1 • J* t I 
M3I*A( 1*1 » J| 

60 TO 30 ' , „ 

30 * locate intersections 

" ?? ICONvIi.LT.AMlWl lR(M| .1 1 M* I I I . 0 R.C 0 NVAl. 6 T . AHAX H il Ml •• (H*ll ) I 

• 60 TO 35 
NP-NP*1 

Oft* ALOG I 0 IB | H* » H . ALOfi 1 0 « 8 I N II 
If ( AOS (RB I .6T.ZFRQ) r.O TO 33 
0-0*5 

60 TO JS . 

J 3 ^ ONJ^OOIO.A.HM 1/88 


TX-OSII .Mi *65 l-l.H I 
TV-05I2,M)*C5(2,M| 


. H I *8 
* M I *P 


X (NT* I »NP I -OT I 1 *N| *CT 1 1. I iM I •Ti*f T I I *2 «KI*T V ♦ j I 
Y I NT* 1 *NP| -OT |2 (X |*CT |«» I »X I* T «*f T 12 , 2 »K » *T T *J J 


35 CONTINUE , 

IF INC.LE.ll 60 Tfl 
NT- NT* I 
„ NV«NT»-NP 
SO CONTINUL 
53 CONTINUE 

* SCALE POINTS 

IF JNT.EO.O) 60 To 
AO DO AS K-l.NT 
NH-NV1KI 
DO AS L-I.NM 
MX,L»*XU.L»*SCA( i 
, VIX,Ll-VU,Ll*SCA| Y 
AS CONTINUE 

• PLOT CONTOUR 





*!CftVN**AStil>tT*re*« 

1 T«io*r,Nr»i» 


* 

! 


Si 


it 


100 


2&SUS 1 1 *B?MiU>oCii otoii^ u*» * o.i I ' 

COMMON /TN|/ Al*NI t0UNI»C|4NI»0U«|| t VU«l 
01 MINS I ON I^HNAIAtl iiff Alit} 

ocT*ur> a fi(jri 

iAMMANrjMOMrl/NrTAl.ifl 

o?!!! 8 ^INCf I tV . 

iffillAWfflfl I/Jctaiii 


•» 

>• 


V (L ) ■IAhISa I L > 

00 200 N* I (LAST 

*»° vT|J"«AMMA( I l«C(| | •¥ f f * I I/Of YAt 1 1 


r,NO 



»|CRON*SnS?l I I. SUB I OfJ 


SUBROUTINE SUmION( IHa«I,JMa*I .* .» .VOIST. fSTOp) 

c •••••• ••••••••••••••••••••••••••• •»»••«••••••••••••< 

C SUBROUTINE FOR 6CNrRAT»N6 ION fNPLANT OaVA 
C •••••••••••••••••••••••••••••••••#•• •••••••••••••••! 

IMPLICIT 00U0LE P-CCisION! A*H«0»n 
COMMON /C0N/6IKIU)|.|M) 

0 IMCNS20N VIII 
ISMO 

jHa JM A I | fZ 

irlYljHItlToliPil) 1 Si* I 

ico ms.miiK«rts:«'*.' N "- T 

!!i»T.m!..s"»»i»r*c«*iyotTjc«T6f .. ... . 


• MORON IMPLANT «V'cF|Oc|»* KEV* ,»| . »CSTnP • *,C|S*A) 

111111:13:1! SS«5SN?;!„ 


MORON 

|P|I55»! 


OO I J-l cJMAXI ... . 

CRI 1 1 cJ>"CMAK«OCXRl-0.»»llVO|ST«y|Jl-RPI/ORP|»**» 
00 2 J"I*JMAX| 

00 3 | • 2 , K 

CRI I I « JlaCMl ( I «<i3 

continue 

RETURN 

END 


M|CR0N*SnS2( It c«BC 


SUBROUTINE *BC J IM*Xl cjMAXI ,X,V,LM.T|M c»K,OF I cTEMP.TMAX.KOTI ,X0, 
• YDIST 1 1 AHfiNT I 
M 

SUBROUTINE WRITES TRANSIENT DATA ON DATA PILE ON UNIT f! • * 

. :: 

DOUBLE PRECISION cR I # X • V • T INC uTENP t TNAX iMfl •X0» VO 1ST 
COMMON /CON/CBliAMcM) 

DIMENSION XCH if 

dimension xi 1 1 1 y 1 1 ) 

DAT* IDOL* IBLNX/| H « r, | M / 

KOT I .KOT I ♦ | 

NARK-IDOL 

IEIKOTI *CT. I \ MARV"IB| NK 
. WRITE! I I .100) MARK 
) rORHATI IHO,*M 

IFIXOTI.ST.il GO |0 2 
WRITE ( I I cl 00) Of I , TH A X * TEMP i Vo I St 
WRITE I I I * 2 0 0 I K * IMA* I . .(MAX I , I AMONT 
I FORMAT II H *6110) . 

WRITE III ,1001 | X 1 1 I . I -2 * | MAI 1 1 
WRITE! I I . ion) lYM 1 ,1.1 . JHiXI ! 

mark-ihlnk 

Write! 1 1 cRoni Mark 

CONTINUE 

WRITE! I I *200) LN 
WRITE! 1 1 * IOC) T | Hr 
00 I J-JMAXI *1,-1 . 

WRlTEIIItlOOl ICR | I I .,1 ) • I • I HA* I I 
WRITE III ,1001 (XOI I > . I »i • I MAX I ) 

FORMAT I 1H .SElS.Vt 
PRINT 300, KOT I 

FORMAT I I HQ , I OX t *Wof I • *,I|D/) 

RETURN 

END 


a £ ^AL page is 

WPQ0 *WALny 
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NICROR* 

1 

2 
J 
r 
% 


Sr»S2< 1 1.« 1 5Ur „ 0UTine mi|HAXl,,|HAXI,* l Y,LM,TIMC,K,OM,TCHP.TNAA,KOfl,»0 

€••••• • Y S ,s !5»5*2*Il#* a •••••••«•«••• ••••••••••••••»••••••••••••*» 

C SUBROUTINE REAO* TmE TRANSIENT OATA MU 0* UN|T || 

C ooUrle precision f») , x*y,tihe ♦TfHP,TNA*,nri*so»YoiS? 

fORNON /CON/CRl I4R »AR| 

01 HENS ION XoT I ) 

n7:l«ri... 0 , . 

_ READ! 1 1 .2001 K » IHaX I , jHAA t , I AHONf 
198 eorhat; ih ,$cis.ri 
200 rORHATl I H iR 1 1 0 ) 


_ READ<li,200t K * IMaA I , jNAXt * I A! 
198 EOHHATJ IN , SCI*. R) 

200 PORHAT j I H .A 1 10) ...... 

RE AO 1 1 1 « 1 00 ) (X(|)il>ailHAII) 
RE AD < 1 1 , 1 00 1 I V 1 1 ) 1 1 ■ 1 »JHAX I ) 
. if AO C 1 1 • 1 00 V MARK 
0 CONTINUE 


. RE AO ( I 1,100) HARK 

0 CONTINUE 
READ! I 1, 200) LN 
NCAOl 11,100) T I ME 
00 I J* JHAA I » I *■ I . 

1 READ 111,100) ICRI | I ,J||t»2,tHA 

2 cSnTINU^ 100 ' ,X0( ' > ' t "** INAXI1 

ROO PORH AT I 1 HO , A A ) , 

APPLY PER I 00 | A B.f. 

IMAA2-IHAAUI 

1 CiS I { mAA 2^ jUcAl ! |H aa \ J) 

RETURN 


IMaAI ) 

M A A 


RETURN 

CNO 




1 ': . „ 


NICNON«SnS?l I | 

A c* 

i c 


* 

7 


3 

14 

ii 

i 

II 
; I 

4 

i! 

i> 

! 

|? 

44 

34 

*! 

31 ! 

» 

47 

44 ; 

•0 

ti 


.PLOT 

5UH00UT | HE ploticx.iha#! .jhaxi ,u* 


SUMOUTJNE PLOTS THjLTBO DIMENSIONAL PROFILE 
IN THE OUTTPUT PHInTOUT IF IPLOT ml 


,mk> M0,,w 


#»ME»lH •lNf*tN » 
( I Hi # lilt 1 1 M I * IMS* f #• I M / 


I 

a 


DOUBLE PHCCISION c»l.rS 
0?HCNS |0 W°STHblI?| , (iUi*INC( 4J> 
oi?J N llHDL)lHi»lH 1 1 HP f I H •iHCilH • 

WIS.iS&i? • ! " y>, “ ' ‘ 

I?i*Lod!o( 2.21)111 . . 

FACT* (ALOfilOUS I *«T )/ (FLOAT IKP*I II 

90 2 J"JMAX | • ( ,» j 
00 I I *2 « | H AX I v, 

IFlCBl < I , J> .LE«0.9> SO TO 3 

« ■ i I alos » o < c, n < 1 . j )> - a T j /r AC T1 ♦ | .4 

CONTINUE. .. 

IFIX.LT. I| ULlNCf * )»OOT. 

IFlK.OC, I . AN0.X.LF*XP1 ULINtl 1 1 ■*>VM§L ( X I 
IFIX.fiT.XPI ULlNfl IlapTAR 

wspiiirt |L> • (IJLINF< M*|a2.IHAll !*<ULfNC<|!tl«tHAI.2 

POINT 12 
POINT 13 

POINT I4*STHBLIXP I iCS 
XOaKP-l 

00 *1 | • X 0 • | , » | 

CBLa( 0. •• 1 1 |w| )*FACT* aT) 
u CBHa 1 0 • ** I I «F AC T ♦ aT I 
4 PNINT IS. SYHBLI It .C bmiCbL 
. POINT IA 

I? Fo 3 hS ? 1 1 Hi | • ! ndfi#* ?? 'rSNCENTOAT t f>N 

• I HO, N III HP ».///) 

FOOH AT I I H 
FOOHATI 

•//I HO ' 

• 1 1 1 1 H 

|4 FORMAT I IM 
»* FORMAT ! IH 
. «o» ,r 10.31 

I* FOOHATI IH , T I 3 , 1 H . • T 2f t * BELOO I.O90EII*) 

OCTUON 


• • 


II 


II 


PROFIl E T / 



END 






NfCR0N*SnS2lti 

2 c * 

i c* 

i c* 

i h 

j p 

i 

ts i: 

t « 

j e* 

if n 

is « 

I! ?• 


27 

it; 

ft, 

30 

1! 

j*; 

»*! 

;? 

35 

3 } 

n 

*• 

«tt 

so 

I 

u 

I 

st 

40 

*11 

it 

*s 

**: 

tr 

tt 

tf 

70 

71 

72 


.CONOEP 

FUNCTION cOnOENH.N,| ,J,M|N.Mi*,rnNVAD 
• LOCATES CONCENTRATION CONTOURS ALONG EITHER 
A VERTICAL OR HORIZONTAL 6N|0 LINE. EITHER 
LINEAR OR L0 «R|THIm|C INVERSE INTERPOLATION 
CAN tC USEO. 

•Cl NR, CONTOURCO tOlNCNSiONCO 
t'N CALLING PROrRAN 


too 


♦ 


* 

* 


J 


CONV AL" 


«?« 

NON-ZERO V«tuc specifics ARID line 

■ WH t « 


HIN, 

Ha* 


• CONOCP- 




I OR J 
CONOUR 

MINIMUM ANn HAKlNUM SUrSCRIPTS OP 
GRID L|NE T? re CONTOURED. IMaV tc 
ZERO tf ENTIRE «RlO LINE IS To IE 

iosi?|nN®fir C0NVAL ON sJtlO LINE. 

VALUES CONftEP 2 (tSrNS^A VALUE^O^ 
ZERO. 


>»« 


# 

t 

♦ 

# 

t 


10 


IMPLICIT OOURLE P* EC \% l ON U*M 9 0»i I 
COMMON /CON/A(AN.A 4 0 
00 100 1 • t t H 

00 100 J*l i N 

If < At I i J) • LT 0*1 00| All »J)*l#D*|00 
| LOG* I 

tr t I .NCOI GO TO 10 

1 INC* I 
I M I N*M IN 
I MA X*M AX - I 

If I MlNoEQ.O) I M | N* | 

rr i NAXtCO.oi ima*»m*( 

If I JfLTtOI I LOG** I 
JINC *0 
JM I N* | ASS t J I 

jmax*jmin 

GO TO 20 
CONTINUE 
JINC* I 
JM I N«M I N 
jMAX*MAX-l 

tr ( MlNoEQ.O I JMf N* | # 
tr ( max.eq.oi jma**n«i 
tr I 1 1 LT • 0 I 1 LOG** I 
I I NC*0 

I M IN* I ABS j I) 

| MAX * I MIN 
CONTINUE 

,IMak 
00 NO JJ«JH|N,JMA* 

200 TO NO . 

If ( ILOO.LT. 01 GO lo JO 

CONOEP* I (COnVAL-AI I I ,jj| l/IAf ||4| |NC.JJ4,||NCl-At|tiJJti| 

1 4 f L 0 A T I I 1*1 |NC>JJ.JINr > 

RETURN 
30 CONTINUE 

f ONt OG" ALOij 1 0 ( CON V AL I . . . „ 

CONOEP”! <CONLOG-OLOGinU( 1 1 *Jjl 1 1 /IDLO.Int At 1 1 ♦ 1 1 NC . J J4 j y NC I I < 
IOLOG | o( A I I I . JJIM |4f tnAfl 1 1 *| I NC4JJ4JINC I 
RETURN 
NO CONTINUE 
NS CONTINUE 
RETURN 
end 


20 
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Cm# 

CM* 

CM 

CM 

CM 

CM 

CM 

Cm 

CIM 

CM 

CM 

Cm 

Cm 

CM 

CM 

CM 

CM 

CM 

CM 

C«t 

CM 

CM 

CM 

CM 

CM 

cm 

CM 

CM 

itM 

Cm 

iC»* 

c#» 

CM 
!C.M 
C •» 

cm 

CM 
Cm 
Cm 
i C • • 

it** 

Cm 
Cm 
|CM 
CM 
CM 
I CM'. 
CM 

jCM 

C** 
CM 
|CM 
c** 
CM 
Cm 
Cm 
Cm 
CM 
C M 
Cm 
Cft 
CM 
CM 
C*« 
C«* 
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mClIH tl •» 

rou S IL ICON ON IaPPHIAI •• 

ON M* M 

M 

MM«M«MMM«MtMtiMMMMMMf«MMMMMMM»M«MM*l«M 

• •••••mmmmmmmM**Mm***#**mm mm#*m*#Mm*####*l* Mm 


solution Of o.rru S .0N i j ! . 8 LfN i ,., 


OaTa IS RrAP fNOM OCCK IN FOLLOWING SESUCNCCt 


FIRST 
LIST . 
IFILE 


LfILE 
IRLOT 
I READ 


JSYEP 
I HA* I 


CARD! 


TO BE SKIPPED »HILr FNlNjINO 

ON FILE ANn ALSO Tfl L 0 CATC 

CONTOUR RO^ATinN.CONfoUR FRON? KOvER|N?_P 
WRITTEN ON UNIT f AND CONCENTRATION 


„ TIME STEFS 
PUT | TO NRilE 


ON 


RRITTEN ... 

IF TIME STE. ... - 

TO R( OT PROFILE PRINT OUt 




AND cONfEN 

UN|t II* , _ 

PS. TO BE * * IRPEp_RH.ll C RBI T INC ON Fli t 


- i or 

-put 

- PUT I TO READ data FROM .FILE. 

- PUT » TO Rt AD | On IMPLANT DATA ANn TO 
DO REDISTRIBUTION. 

. • of oaTa steps to be read from file if ..... 
jhaxi - * or 6Rin foints in * ano t direction 

RESPrfT | V ELY • CHECK DIMENSION BEFORE CHANGING 
‘“ *110 


format fiflo 


s / 5 «““ S6?°i ,r const. source oifp. is oes 
- PUT I I r R f ft 1 5 T R I B U T 1 0 N IS DCSIREO. 

ORINt’-^pIt ! FOR T ipO S fRTSTAL r oRipN?A?ioN n 

: § *88 ! ? mill 88IS8WI88 

AMBNT - PUT I FOR ojY ftXVGEN 

. PUT ?, FOR rTEaM 

- PUT 3 FOR MlTRftGEN 

FIELD JllO 


DESIRED* 


CONCFnTiIatION/I *E|R 


third card: 
csur - substrate 

CS - SURFACE CONC. ... 

temp - tempcraturf in deg. cent. 

tma* - NORMALIZATION T j me in second 

THIS HAS NO EEFfCT ir LaMoa IS SPEflFlCO AS OaTa 
K olST. To I ST - WIDTH and THICKNESSIIN MICRON) OF THE TBO 
DIMENSIONAL region in question 
OKTHIC - RIDTHIIN MICRON) or the OKIOE in the region 
dflt - normalized time STEP . 

FORMAT FIfLD - *F 1 0* 3 


lMTYPE C - R S^Ef irV TYPE ftF IMPURITY 8j PUTT 
NO SPEC, i S NfCESSaRV IF |T IS B 


|NS N 
BnRON 
OR ANT 


IMPUTV -PUT BORON. ARSrNlC* PHOSPHOROUS 

other name* 

ACTIVATION ENF RGV OF THE DIFFUSION 

if blank ano rORon d i ff * . Data is supplied 
DIFFUSIVITT const, of the IMPUR ITT • 

IF BLANK aND bORON.OIFF., OaTa IS SUPPLIED 
FIEl 0 AH.rAH.PF 10.3 


Ea 

01 


ON P 

INTERNALLY 

INTERNALLY 


ON COhHENT TO BE PRINTED ON TOP OF 


AT NH,cH s SlMUL fi TION 


.... STOPS 
ANO^SafFHIRE 
PREDCp» 


FIFTH CARO • 

10 - lOENTlF |CaT .... 

PROFILE PPlNT OUT 

ITEST - PUT o TO RrAO'LAHDA FROM DATA OCCr 
CSTOP - CONCENTRATION/! .EIS AT JHIC 

WHEN THE lfFT fNO CORNER OF 

INTERFACE REACHES THIS ValUEIOURInG 
FORMAT FIfLD - I SAG , 1 S.F I S.f 

SIXTH CAROlUse IF .iSTp*0)5 

ROTEMP - REDISTRIBUTION TEMRERaIURF 

RoTMa* - REDISTRIBUTION NORM. T IMF 

REDISTRIBUTION .FINAL TImf is I* 

HOOLT - redistribution norm, time step. „ ... , 

XOA - RCOlSTRlBUTlftN INITIAL OXlOF TN I cKNrSS IN CM 


• R 
RR 
RR 

RR 

RR 

RR 

RR 

RR 

RR 

RR 

R* 

M 

M 

M 

M 

M 

M 

It 

• • 
M 
M 
M 

M 

• • 

M 
It 
• • 
• • 
M 
M 
M 
• • 
• • 
• * 
M 
• • 
M 

M 
• • 

■ M 
M 
It 
M 

•!» 
*]• 
• * 


•o’ 

jl 

it 

1 

i 

I? 

!!! 

III 

!|j 

10 * 

t if 

ii 

it 

117 
III 
1 19 

ill 

H 

lit 

»• 

! i 

|| 
A 33 

iji 

I 34 
1 37 
1 31 
1 3f 

• 2 ° 

{I 

|H3 
f 11 

I IS 
1 14 
1 17 

i;5 

18 ? 

|SH 

is* 

t« 
1 $♦ 


•herf sukf, cone. mas cs 

Xon REOISTRIrUT | 0*1 I NjT IAL OK I or TMIf KNr*S IN CM 

•Mctt mint, .has the thick al i Dr _ 

CH - SEGREGATION COErflClMT* •ENMaTKO |NTCRNALLT 
IMPURITY IS rORon AND NO VaIIIC It t|V(N 
FIELD 4FI0.3 

seventh caSdiuse »hen itest-oii 
u "“* - r,,., 

f AA M jBp£*li?iV| 6 N 'Eij !« ' * 

RR*5 RaNGk'oe DISTRIBUTION. HEAN VaLUC (IN MICRON | 
ORP • STRAGGLE. STANDARD DEVIATION. (IN MICRON I 
VOLT - IMPLANTATION ENERGY LEVEL IN REV 
MELD 3ClS.4.plO.S 


If 


•i 
• • 

• • 
M 
• • 

• • 
II 
M 
II 
• I 
II 
II 

II 


, S? L,c i T «S?H 9 L ?R- p,l f C| * 1 0N • io-i i 
INTEGEK OKINTtAHAN* 

DOUBLE PRECISION kSiNfiLAMOA 


C 

c 

c 


c 

c 

c 

l 

c 

c 

c 

c 


Citiiif i*<H » 4 CAM f *4 1 • 1 01 i 1 1 


DOUBLE rWECISIUN KO»Nf ILAPOA 

common nn\i au<d •oiahj icui) iVIM) 

§IHENSI0N°l(l4N?! < $(ii| . CBIAl.All, CltUNt. . 

DIMENSION IHTyPEIll. fNPUTVm. IRRONU). 

DIMENSION ROLASTlAll. R0I4RI, DRRIO lAII.VPIAII . 

DIMENSION |NrT(7I«I0rNT( 1 1. 1 OOIItl.INITltl.lt TH III 

orr inc computing eiinctiOns 

?8?5i??5!8^??2?Ui?t5i8SSf!:ii?VR00T 

CALL CRTRAN (T.IDaTE.iTIMEI 

read simulation iniTiai data 

READ 140, LIST.IElLEuFlLC.IPLOf*, |REAO, JsTEP.|MA*| .JMaRI 
RE AO 340, JSTP,OR|NT,*MdnT 

iEii*St8: c ?s8i.«?^.T»» •XDIST.VDIST.dRTHIC.OELT 


READ 320, 
READ 3B0, 


!MTVPE.|MPutT,EA,D| 
ID. tTEStUSTOP 


IE (JSTP.GT.6l Rf»0 I S. 370.CRR*3 1 0 1 BdTEmP . RdTHax , ROOL T , lOA . X abTcm 
IE < ITEST .EO.O) READ |$,370tERN«i)l0t LAMDA 

define data: 

RP - HOLT 2M ANN • S CoNSTaNT 

Jt. 1M ■ I Of aLLOWarLE. ITERATIONS 

DVLIM - ALLOWABLE rgROR IN CONVERGENCE . __ 

TimAK - NORMAL I Ef.o Time at rhich SIMULATION stops 
Ear and DIB - BORON OIEFUSIVITY Data 

OATA FAB.OlR /3.SRD0, 1* 1700/ 

data ibron AlHBORO.lHN .ih ,ih A 

data nyp.ipt /inn. I HP a 

OATA KB AB.AIAD-G/ 

data jlim.dvlim /aO.i.d-6/ ,tihax 

data IDOX/4HORT 0«i4hyG£N a 
DATA INITAAhNITROG.AHeN a 
data ISTMAAhSTEaM a 


A 1.00 A 


,MTYPe-IPT 


Vf 

IF 
IE 
• F 

IF 


I ER» I 


ICR»i 

!£*■( 


IE i I lHPUTY.ro. I bRoN) 

|ER«0 

I JSTP.EO. I.ANn. IK faO.EO.OI 

I I MPUT Y . NE • Ibron • AMO.CA.CO.0.00) 

( lMPUTY.NCt IRR0N.AN0.0I .EO.O.DQI 
i lFR.EO.il PRINT 1R0 
( ICR.EO. I I GO TO 3 | 0 
{ I MTYPF.rO .1 P T V IT Y £ E ■ 0 
I lMTYPE.ro . NYP . ITVPE- 1 
( | MPUT Y. cO. IBRnN. ANO.EA.EOtO.nOI E A*E At 
I IMPUTy.ro. 3BRftN.AND.Dl .co.o.oni 0 1 *0 ! • 
PRINT RIO. ID .. 

PRINT 330. IMPUTY. IMTyFE 
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ill 

58 ! 

*07 

*01 

*0? 

*10 

III 

*13 

if! 

fit 

IS! 

in 

iff 

Hi 

iff 

*** 

if? 

if! 

if! 

lit 

fit 


C 

c 

c 


GO TO % 
0 ei.K-A. 0 -H 


c 

c 

F: 


NAMELIST /PUT/Ll 5 T,|rlLE»LfSLE« I PLOT .IHfAD.JSTEP.IHaxI . ,|HAX I « JrTP 

5 :?!} 8 J:cXI:Sm! " " ' 

•HITE 1 4 • PUT I 

OtroiNATfON or cohrUt imo« geometric AMO PHYSICAL r«OAMfTc«f 
Iff JSTP. EG. 0. OR. jRTP.rtT.lt ICONO-I 

iriAHONT.NE.il 90 |ST«v0|ST«I.0-H 

ISlISSaiJfEr 1 ’ 

COT 1-0 
ih ax«i max t -i 

| M A X 2a I MAX I . I 
jci iHAXi.sT.il 

|(VaHRN?. NE.i) 

60 TO A 

* 0eLX«XD|5T/rl0AT(|HAX*-3| 

* 0CLY-YDIST/FL0AT(,|MAX) 

CSU0*CSU8*I ,0|5 

cSTor.csTor.i.oio 

LAMDA-LAMDA-I.O-i 

es-csAUOia 

T-TEMP^271. 

RINOOHlXOIST-OXTHtCUl *0-5 
GHaCM 
r.H» i .o 
THX-TMAX 

iri AHBNT.EO.i l I HrtTt I |-ID0X( I | 
triAHRNT.EO.I I IHaT(2»«|00XI*I 

lfuS»NT T : 5 »: 5 i !tt?Ti?:r« m „ 

!SS5«!TA8l T '** 

lMORlMT.ro. 2 1 I OCNTa ♦ j 1 0* 
iriOriNT.EO.il lORNta* III* 

calculate distance x ano v 


io 

*o 

iO 

21 


00 10 laZ.lHAXl 

Xltl-I l-2»*DELX ... 
IF (XIII .LE.HlNnO) 
00 20 JbI.JHAXI 
-I »»DE 


If • I 


Y4JI-I J-l I AOEwY 
00 30 I *2 . 1 HA X I 
XOL AST ( I I aXOA . 

00 31 JbJMAX I .1,-1 „ 

YfiJlal JHAXI-J)*I I ./n.OATI JHAKI I 

srecirY predep cowriitinNs 

00 NO I ■ I . IHAX2 
DO NO J-l « JHA X I 
Cfll I • Jl-0.0 
CHI 1 I .Jl-0.0 
NO CR2I I , Jl-0.0 

00 SO I-2.K 

CRt t • JMAXf >bCS 
CBI ( I . JHAXl l-CS 
SO CB2( I • JHAX I l-CS 
TIHE-0.0 
lHaO 

specify initial proFIlf 

IFI IREA0.EQ.2I CALL SUB I ON I I HaX | .JHAX l » K ,f» YOfST.eSTOP I _ 

IF IIREAO.EQ.il C»LL xYZ < | MF I , JMF 1 . X , Y ,L M , T I HE , K A .OF A • TmP *T JHX , JR 
ITEP.XOLAST . YOIST.aHBNT* 

IFllREAD.EQ.il T ImAXbT I HE ♦ I .0 
STIHEbTIME 

IF IJSTP.EQ.il ICONP-2 

NEAO CONTOUR FRONT HOvrHENT FROH RAYA FlLr IF |REaO b j 
IF (IREAD.EQ.il C’aLL fRONT 1 1 H I . ,|N I ,L A *Kr . THPR , OFB .OUT ,LM . I .COT f , « 
' IF S » IrEAO.EQ. I I PRINT HiO. LH.TlMC.TNP.OFA.TlHX.lMri.jHrt ,CA 
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**0 

\l\ 

S 

IV 

II 

IS 

III 

B 

III 

ii 

us 

m 

III 

ill 

ill 

ill 

jBt 

*21 
*2f 
*2* 
*»*• 
*» 5 

111 

IP 

SOI 

SOS 

SOB 

SOS 

304 

307 

SOB 

309 

510 

511 
S I 2 

ill 

SIS 

lit 

ii; 


i 

e 


i 

C 


c 

c 

c 

c 

c 

£ 


C 

C 

C 


i r p*!nt ‘♦ho* 

If I l*^t • f A*K?«oi?,|M? I .NE.jHAli.ORaKA.NE.K I I (Ml 
IF ( IMl.Nt. IMAM .ON, JMl.NC. JMAMI .OR.KC.Nr .K.OR.L A.Nc.CNi |ER«| 
if urfttEQtii mint «tsO 
If (lEft.EO.I | 60 to 3(0 

ST40T SIMULATION STEF 

so |g r ( tgONO«CO« I t 60 To 70 

te&mw 

tma kIS8?^ax 

If I AMMNT .NE.St Oft T*OfLT«TNAX 
If ( AMRNT.NE.3i TIMAX«tHaX 
T I ME>9» . 

70 CONTINUE 

Calculate ni 

T|MAX>T IMAX-OCLT 
CALL PAN AM (NI iT) 

Of>DI>OEXP(*EA/(Ki»*T| a _ 

If (ITEST.NC.Oi nrI»Or*TMAX/l |t0IST9l.0-Nl**|| 

I F ( AMrNT * NE « 3 ) OFt*DF» u 
If ( I TEST *EQ*0 I Or I “LaMDA 

If i I CONO «EQ • I I POINT SBO # ■ 

If (ICOND.EQ.2) PRINT *50 
PRINT 390 « Nlj 

PRINT BIO. 


lfU8?y ; 


TEMf.Dfl.TMRX 


P* lor |*OELT l/|0ELT a *2l 
• P(On»OELT)/lOEL***2» 

PRINT INITIAL RROf |LE 

- c AVV- oU ICV , 2-, * 5 aX * 1 • jHAN * *K rLPll JJ I Hr • yd 1ST • |f»R I T | HC 1 « Ok AST * 

• FT IHf ft OT J HC • AfVoMf I 

If < iRCAO.EO.I.ANO.IflLE.EO.l » CALL, ABC ( I MAX I . JMA* I • I . Y .L* . Y | Mr »* 
I ,0f| .TEMP.TMAX.XOTl .X0LAST,Y0 |ST,aMBNT» 

*M» 1.0 
CCal.O 
RR» I • 0 

IF (ICONp.NE.il CALL oXOATa I AHrNT.OR INT.T.BR. 

JCCtAHoKB) f > 

If I ICOND.EQ.2. AND. OM.fO.O. I CM MM 
If ( ICOND«E0.2*ANO.fiM.NC.O.) CM a«M 
|f( IC0N0.EM.2I PRINT B|| . I HOT . | ORNT . Br .CC .CM 

START TIME STEP LOOP 

BO LM«LM»I 

STORE N TH. RESULT FOR N.H.S. .WtLI NOT BE CHANGED ^UR|N6 I TER. 

lf7!"oNO?EQ^) ST|ME»TlME 

?f 1 ICONo5Eo!2 , ?%T|Nt«TMX 
otime«RoTmax«(T imc-stiHe » 

If IAMRNT.NE.-S) OTjME.tIME 
00 90 I ■ I • I MAX 2 

00 90 J>t . JM A X I i 

CB2 t I . J > -CB | f | • J » , 

If IICONO.EU.il no To 1 10 

CALCULATE OXtOE THlCRNrSS 

|F( AMRNT.NE.3i GO TO fS 
, 00 9B I >2.1 MAX I 
9** XO ( | ) ■ XOLAST ( H 

CO TO Ito 
95 PS>BB*DCLT 

OS-BR/CC .. „ 

00 100 1*2* I HA A I „ 

„ XO 1 1 l"XOL AST ( I l«PS/|2,»X0LAST| I UOSI 
100 ORRXOI I )>|X0( t I.XO laST(I) l/DELT 
AO « I I •KO I 


90 


181 


lii 

• w 




tf<<YOIST-e.*l»*»Oll»»,U.O.Il> HI «0 TO 1(0 
1 10 CONTINUE 

OeLT.(| | . / PLOA T UN«XI|.lTOMT.O.NN«ie|.AtT|tM 
truH0NT.Nt.li p.,of i.oelti/|oclv**2i 

JM»0 

ST ANT ITERATION LOftP 

I to JJ* 




C 

it 

c 


ir I JJ.fiT e JLIH) POINT #00* CTf 0 
If I JJ.6T . JLIM) Oft TO IfO 

TOANSrcN #OLN* VEctOO rNOM LAST |T(*A* PON CAL* Of « 

00 ISO .l-t»IHAXt , 

. .A 00 l?0 U-I.JHAII ... * 

no col i .jueoii t ,ji 

cal*# 

00 1 00 1*1 * I MAX 2 . 

00 INO Jfl • jhaxi . 

OA-OATIOICM I ,j) .fSUi.Nl » 

RO"ROOT | 0 A I 

(NO 6(1 .JUFUIRa.rO) 

SOLVE SYSTEM IN V nlRCfTlON 

00 220 I • 2 * I HA* I 

Calculate the coefp* in y direction 
sr»o. , 

00 ISO J"? . JHAX 

6IM6I 1 , j|*S| |H ,.|) 1/2.0 
G2MGI I i JI*6I ||J«| 1 1/2.0 
63" ( 6 ( I , JU6U-I , j» 1/2,0 
GNMGI I . J)*GI | *J-|j 1/2,0 

1 r ( AHHNT.EQ.il GO TO |S| 

... SMO.NStORPxOI I | • | Y P I j I ■ I * 0 1 • 0 C L T / 0 C L T 
151 CONTINUE 

A I J I ■•P*6X 

M JIM.*e*(GI+63|*P*(G2«GM+SR 
„ _ 0 ( J I ■ *P*G2"SR», „ 

ISO C ( J ) «RHS I Q « 6 I ,#l,rO(l-t . JI.CM2I I.Jl,rO|l«l»U| I 

put boundary condition on y axis 

M7)>A(2U0(2| 

A ( 2) ■0.0 

IF (ICOND.NE. II 60 jO |7q 
I f ( I .GT.K ) GO TO- 1*0 
E| JHAXUE( JHAXI-0I JHAXUcS 

01 jmaxuo.o 
60 TO leo 

140 CONTINUE 

Bf JHAX )»B( JHAX 1*0 | J*AX) 

OtJMAXUO.O 

170 


60 TO | 80 

IF (AHBNT.EO.3) ORPkOiII.O. 
Ma*"DELT* I I . /CH"0* Ns I •DnNXOt l| 


c 

C 

c 


H *N*|»Ha}|HB*6( IU.iHAAI I I 
WK2*HA/<HS*G( I ,,|MAXt » 

§ { J M A K » -0* JMAK * * * 1 1 * * MK 2 ,/(i * *HF 1 1 » • 0 ( ,|NA X I 
(00 CONTINUE 

CALL TNIOaO ( 2,.|HaX| 

CONVERT MATH | X SOUlT | O n 

00 |90 J-2, JhAX 
CM ( I . J I "V ( J I 
I 90 CONTINUE 

PUT NOUNOANY VaLUFN IN Y 

cm 1 1 * I i*cm 1 1 


no I 2S. AI ' PAGE 
0P POOR QUALITY 


182 


If iicono.Nc.n so fo a in 

If (l«fiT*i(l co to ino 

cm 1 1 ijnari i»cs 

tpo S8,I?,5^„.cmu.„HM, 

60 TO 220 

110 CONTINUE v 


; MT OOUNOAfV V*i,Ue4 IN AXIS X (PERIODIC OOUNDAfTl 

. ISO CKI.I IHAX2 •«M , CR| I IHfXtJl 
: CHECK FOR conveRof.nce 

j jF^J.EO.II 60 TO I2r» 

00 I?0 I ■ 1 1 I H AX I 

00 |f^«C»Tl!ij|!(,E.O.O| «o TO (NO 
CTCfMDHAX|(CTfST,r.TCm 

(NO CONtSnUc „ / 

If < ICK«Ni:*0) 60 TO |*0 
JM»JM*I 

do am i*i » i{JaXi 

>«> w 1 ,i5iw:iii«:*.»»i 

; IFtJM.E8.il 60 TO 120 

I MINT RESULTS 

|S»LM/UST*LIST-lm , . 

| If (IS.EO.O) 60 To 2S« 

If ITI)1C.6 E.T|MA*,ANo.CS.Ne»0.0I so TO 2*0 

If <Cml2*l).6E»£ST0f ..ANPtlcONOtEe.il 60 TO 2fO 

If ICRH2* I I •6E*CST0P. AN0iCS*EO*O.O*ANO*JRtf •CO*H SO TO fSO. 

60 TO 2 AC , „ , 

(SO # CALL OUTfyT^I^TiiHAXIfJHAri.KiLH.JJ.TIHf.TOlSTiinrlTlHC.XO, 

CALL sArE JN* IcSUr* r 0 1 sT e JM aX I • X * T |HC tOCLT *0(LT t TcMf « THAI » IflLf. |TV 
I PE • AMRNT . XO ( 2 1 1 

iwiwwiiii.jtii.SHM, 


cmt t ,j.»eti<i*j) 


•*» 


STORE TRANSIENT DATA 


(AC |f IIFlLE.EO.OI 60 To 270 
| T»LM/Lf I LE«Lf ILc-^M 
|F< tT.EO.O) 60 TO 2A& 

*“ I I!fef “ t II AA , :i S|»«S T I ; S »: ♦fef*T I SXAht » 


,*N0. ICONO.EO. 1 1 60 TO (AS 


Locate contour position ano store oaTa 

2AA CONTINUE 

„ CALL fRONT I | MAX I , JMa* I V tH .Kilt HR .Of I .OEl T.JSTEP ,0 *K0T I , KOIST | 
(70 CONTINUE 

60 to next Time stirP 

IF (CRM 2 . 1 I *fiE.csTOp, AND* ICONO.EO* 1 1 60 TO (00 
If IT|Mc*CE*TIMAXt 60 TO (fO 

00 , 2R0*|iIllMAX| 

l#0 oo*?k*8J * ' 1 \ 

(TO IF I J’STP • 6T • 1 • ANf( • I COND« CQ • 1 1 60 TO SOO 
GO TO 310 

300 timax»time*m 


■ ! V;: 

I . 
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•or 


110 


irUMNNTiNC.lt 
ICONOaf 
fiO TO NO 


TlHANaTlNCNNOTNAx 


Hit" 


.CfiiOt MINT ASo 


||8 { 8 sa» itMUiftibUm 

!!§ ;s: 5 i; !HWMs« «««««.■ 


innun m Type # •,*•» 


cyctc •••«»„ 
erctr ••••» 


138 iHAy'.U.riiiti 

»,WJ!^!'!!MII?inin«^ until! . 

L8 N iflS P Cr^LFflaiS2l / *iS°^i ,, Al32 , 5e. t . , iK? «0lUT|0N»7/lH •|0»»|i*R//l 

**V? StJi«?{ V S?NS!l{? 5 ;%i«!{ 

f I// 1 N *» 0 ?.!?aH»a*: *r* 22 **?£^-* ;unt 


Nil 


•//IN .IOX.*LAMOA," LtNOA**> • Ori«TNAX/| VNAK«I •C«A**»M//I 
I'UHOA • •,0l0tS,»X|tr?R NORMALIZATION T|NC a ' »C 1 2 . A . Zl , *1 

s y(i «;*!'» !«■ «»*■«« 

f»o»iiiXi*5_?-!i?tS«i./_ _ . . 


flO rONNATIIwp.lO*. |N| Tl At 


,CC.*| 

ONICNTITION • 




|j « islf ^5 ^ i ?is !y^AAis«*» 


oriLC At UHC STCr a 
*■ •Ntt/j 


• • .010.3, 2X/IH 

MoVoRMAT 
rue 


.laX.tNONNi T INC a * 


i f 


■Mi;!?! 

I I 0 . 2 H i • JM 


* $0| Ot 3 *XA . * INAX I 

. . , , . ARC OlTA (NCO fROM THr OATa . 

f ••»l//|M_i|f|Xi*T(NP a t,fiOil/|N 1 10* • *LAM0A*«2 • *»DtO. 

• * 1 10/ IN .|0«*MNaM a ••tt/lH .tftXt'UMAtT 

NiO fONHAT - ( I HO. I OX • ••••OATA INPUT fRROR •••*/IMOi*«.*RUN«TNtAN TcRnIn 
lATCD* | * * "* 


il/tN ..OX..TINC NTCP 


(NO 


•PRT *0*» tSMRC JN • .HA | N-PtOT 



C f poor OH?* is 


«« 0 IY 








A|CA0N*5fl*3!,l!.SHNrJN . . . y 

SUIROUTINC SHREJN*CSUa.VDIST,JHaX| ,K,TIHr»OElT»OCI.V, 

• TCMR.TMAX.irlLC, I f VlE • I AMINT • XO | 

(•••••••••a «••••••••••••••••••••••••••••••• ••••#••••••••*•••#*••••••••• 

C SUBROUTINE for CALCULATING SHEET RESISTANCE ANO JUNCTION 

INFLICtT OOUSLE RalC ISlONI I-H?0*n " * * 

COHNON /CON/CIM*a 1 » k 


l 

I 

s 

A 

f 

0 

•I 

It 

H 

4 

; 

• 

It 

R 

23 

It 

{ft 

24 

it 

2§ 

It 

II 

11 

37 

SI 

39 

3° 

5> 

*11 

AS 

HH 

3* 

AA 

AS 

AA 

A? 

SO 

it 

ii 

it 

ftft 

if 

t? 

AS 

AS 

At 

AS 

AA 

ti 

n 


IF ( | ANSNT (NC*S I 8sfL»©rl/YO|ST 


V«J«0.0_ 

VRJIaO.O 

00 A IJI • UMAX! 

ij-jhaxi-Ni 

! ! amrnt Inc . s i ’ v t i i ■ t i * i » • it • /float ijnai I • 1 1 * 

jN|a| 

|CKa(l 

DO S jaJNAKlili-l 
I aJMAX I *1 »J 


IficA!?! IIC tIcSUR •ano.ick.co.oi 


C 

c 

c 


... jNI"| 

iriciii u I .gt.csuri ic*"t 
IF I cat 1 1 i .gt.csuri JNal 

IF? Jill »eS. I loA • J* »?8» ( JN AX I •! I .OR. jN.Sli. jNAX| | SO TO V 
VO J» I VO I ST / UMAX I - lift (CO 1 1 JN I •CSIlO I/ICO 1 1 JNl *C0 1 <JN*| 1 1 
IFMAHANT.NC.3I V«J*ViJ/YOlST . . ■ 

VOJt»OCL*(CSU»«CRl (JN| > l/'CAII JNl«l >-C0l|JNt I I 
K«0 

f ICHAaO.OOO 
0*0. DO 

START INTEGRATION 

(•Ml 

IFIK.6E.JHAX! I GO JO 9 . 

Q*OA ( CB 1 1 K*| ) *CA| | K| I •OCL/X.O . 

|F<K.LE.JNi.CR.K.«EijN> GO TO I 


1 v I k ! I • us iff ti is i i/a. ooo i 


# « « W T i f nnv 9 W If » f 1 n ^ 

1 1 62*6 I DABSlCNET (If * 

A I GH A >S | GM AM V I Km I 
GO TO I 
CONTINUE . . , . _ 

IF< JN.EO.JHaXI) VmJ-O.?. ^ 

irl JNl .EQ.I .AND.C nU JN j I .GT.CSUll VRJ 1*0.0 
SI6 A*g(DARS(CNCT | ,|N* | I » • | TYPE » 


9 • w p n •• • w n w « i w it i i *r/ : • w r * • » w r . 

RIGRaGlOABSKNETuNI I MlTVRCl 
$ I GCaG ( 0 ABS ( CNfET I ,|N( ♦ | I i , ITVPCl 


VJUNCaVl JNUVAJ 


Y*J|aV< JNl IMRJI 

5isr:53J5S| r :Si‘ ,J - . „ 

>0> ?8«Ut> 5E S 1 • jSS/ILS^'TS •!?*•• 


VRJI|*(SISAaSIGC!I 


•sx. 


5 .... .fois.f • * 

!.0lS.9j3X. 


CM* (SX./. 


SOI 


. .‘SMEcf NKSI5TAWCC ■’•OIA.t.5 
/SX. * INTEGRATED I HFUR I T V ■ * • 0 1 S . 1 1 
•R I TC < A . 20 1 I JN|,jN . 
rORHAT(*0*.RX.*JNi*.2l*/> 

AFACvl .0 

IFUAHRNT.NC.3I Ar ACa f VO I ST-O. A5«X0 1 / 1 .0-A 
VJUNCaT JUNC.AC AC 
RSaRS/f Af ACM .O-Al 
QaQMP AC* I « 0*9 

tcruiM 

CNO 




MI(»OM«SnSf 1 1 1 .MlN-PLOT 


I 


ft 

t 

li 

It 


[ft 

Ift 

II 

It 


n 

s 

I 

if 

II 

IV 

ll 

3 

;* 

tft 

ft e 
:> 
li 
ti 

ft? 

»• 

« 

*i 

a 

t& 

** 

*? 

«• 

*t 

?o 

? i 
72 

it 

n 

?? 

?• 

?t 


•• iSOrONCrNTUftTlON PlflT PAOGRaH «• 

*1*11 clt'll it! I'J"oJ‘{fH?‘5?C*, TO u , 

ttCONO '»*., I IT ?S 


* 







a 

Cft 


OATA IDOL , IftLNK/ |H>( |W / 
DATA M6T/0»0ft?ft/ 
itftO 15*3) LL 
•tA0tft»3) Ilf 
M 3 ) 


FOAHAT 

oo*fo. 

read i? Ml 

FOAHAT I I HO i Aft ) 


A 7 ft, f ItC 

Akkic 


Ti 


tt 


c# 

Cft 

Cft 


aq 


IFMARK.NC* IDOL) flO To ft 
RCAOM. I ) DFl *fNAVtTCMPtVD|Sf 
AC AO (7 *2 ) K, IHAXl , JH anI • IAHINT 
IK* |HAX I 
JK.JHAX I 

FOAHATUH iftClS.tl 
FORMAT! I M fftllpr . . 

ACAD 1 7.1) < XC ( I ) * i*2, f K) 

AC AO ( 7 *11 MC(I ),»•». Jft) 

AC AO I ? » R) HA AX * 

CONTINUE 

READ I 7 » 2 j LH 

RCAOI7.I) TMC 

00 ft J. JK . I *-l 

•*£*?* 7 f * i Cftll.Jlftl.FMK) 

no 10 l>2i IK 

ificai i.Jl.LC.l.r.aoi CiM.ji-i.r-30 

COn T [ NvC # m 

ACAD 17(1) IXP(I).MIaiK) 


|S*L/f IT*I |T*L 
IF 1 1 ft* CO *0 1 GO TO 


II 


if(l,«eq*i .OA.i.ro.cu so to ii 

it cSnfi»4f° 

• INITIALIZE PLOT. 

'.‘tisWlll'WiiMI. 

CALL PLOT! 0*0* THAI *2 »• I ) 

CALL PLOT I XHftX » VHaK »2. ■! ) 

CALL PLOTIXHAX .0.0*2, • I 
, CALL PLOT ( 0 .0 ,0.0 , 2 ,- | I 

• ORA* OX I DC HASH 

IF* Xft-lXf I IK)-2.C-R|/lftC-R 

caLL^plotioxmHax ;a» 

fitt ft 8 ll 8 Sj»«ilU;i:;!iii 

CALL PL0T(XHftX.TH»ft.2,-l) 

X.XHAX 

PAT.ft.-OX 

WWVJWc.' 2 ' 

X-X-0.2 

CALL PLOT X,VHAX*ft;2,»'-|) 

„ CALL PLOT I X , YHAX , J I 
10 CONTINUE 

CALL PLOT 10. 0.0* 0,3 I „ 

• HORIZONTAL grid linfS 

00 20 I - I . t 

t*o,r*floatm I 

1 NAT .ft 3 

IF I I.CS.bl LNPT.*I 
CALL PLOT (C.O.T ,31 

• VERTICAL 6AI0 |.?NCS 
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•e 


! 

i! 

I 

i! 

iff 

|02 

103 

{8* 

}5» 

18 * 

I 0 

Hi 

Ip 


: 

‘ii 

i‘* 

li! 

15 ! 

32 
J 

i»S 

I 24 


I! 


|3t 
I HO 
HI 
1*1 
I.H3 


00 30 1*1* It 

«.o ?2:5 u°mn 
ir », Hooii.ii. fO*o» lnpT**,i 


«♦ 


*.0.0,11 

30 CONTINUE 

t MONJIONTAL A* II NUHftCNS 
[RMP/20, 
f»2 1 


Ct 


200 

100 


V . ¥ 

ill 

ilT 8 |s:t*: 8 ni list 1 


!<• 

B(l;l*/20« 

31 CONtInUE 

m “MiMSt 


nc«MP{io» 
•I) TO 3 1 


laiwtftmilrii. 

V—0.2 


»*M«7*0.0*I 


^o ISfcVi2ot BKMR ’ v * MRt|¥A ^ e * 9 * ,#l 

• VERTICAL, AXIS NUMBER? . 
l,r ( I AMBNT.ur *3) y,k*yc< 
QEYfVAK/IO* 


00 fto I* 1, 1 1 
¥At*Dl Y.rlOaY 1 I 1 


II 


■••3.0*M«T-0.1 . . 

V*YMA)l-O.N*rLQATI f Jl l.0«**W6T 
§c CALL NUHB£*U,Y*MbT,VaL,0.0*I> 

CALL SYMBOL (**.2 ,3.7 *0.1,*- l.0E2n**0,0*B| 

11 sipiii ill iii ill mm 

CALL S Y H 0 0 L I 0 , 2 . H M 1 * o, },♦• 1.QE8 i ♦ « 0 *6*1 1 

Mil I S88 ! !8;f: :8: !f:8!8::i! 

CALL SVHR0LlN«e»3,3*Q, ,2,0.0, M I 

lltt s ti : :2: :S: ; :S:8: : n 
lit! !JS88tlJ:8:{:3;8: 

(HI SYMB0LI0.,S,N..2,*»I ,0.*-l > 

EStt 5!S8flii5!!u;}!«:l:58i:8::i! 

CALL *?NNOL 10.0*3, JiO.*. ‘TEMPERATURE 
CALL NUMBER 12,4,5*1 *0. 2 , TEMP *0*0 ,0> 

CALI SYMBOL 1 0*0*0, ■ *0.1* *T INC SYrp 

CAL * ^ NU M B E R I f » A , 0 • f * 0 . f . T L M • 0 1 0 « • 1 1 
CALL SYMBOL 10. 0,A,B.0.2 # «T|HC 
CALL NUMBER 1 2.4 ,*.1,0. 2, TIME, 0*0. It 
Afi 200 II* 1. A 

CONVaL* I O* 0*«FLOAt 120*1 1 ♦ 1 1 

CA^rBLTCOMICB.CO^.n.U.^, 

CALL PLOT 1 0.0 » 0 • «* ttf » 

CONTINUE 
STOP 
FNO 


• *.0.0,101 

• *.0.Q,J«) 

• *,0.0.1*11 

• * ,0,0, 1*1) 


OF IN 




S0inq Pag* 


not num 
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APPENDIX D. 


FINITE-ELEMENT EVALUATION 

In order to Letter understand the applicability of the finite 
element method, it was applied to a one-dimensional linear diffusion 
problem. This simple problem Is one for which familiar results are 
available for comparison and at the same time taxes the finite- 
element method. In its most valid form, the finite element method 
is applicable to variational problems in which a true minimum of an 
energy-related function exists. Such a minimum does not apply for 
the semiconductor problem in which current flow occurs by diffusive 
and conductive mechanisms. It has been proposed that a "weak form", 
the so-called Galerkln method, be applied to such problems. The 
typical semiconductor problem is a non-linear boundary value and 
initial condiction problem of which the linear diffusion problem is 


a very special case, 
u, obeys: 


In the example chosen, 


the diffusion variable 


4 


1 


u(x »t>lx-o “V * 1 

(a) 

(D.l) 

u (x,t)j ■ 0 

x / 1 x*a 

(b) 


u(x,t)| n - 1 , x-0 
1 U 0 , x>0 

(c) ; - 


u -u * f (x, t) - 0 

t XX - 


(D.2) 


The Galerkln formulation of this problem is: 



\ 

} 



2 


3 


0 1 

t 
k 

Cbl Bicubic (Harm I to) til I f functions 


Flguro D-l. Illustration of hill functions usod 
In f Ini to olamont mat hod. 





where v(x, t) represents a "trial function" which Is used to 
approximate u(x,t). The finite element methods uses a set of 
"hill functions" as illustrated in Figure D.l to construct the 
v(x,t) approximation. Two of the popular hill functions are the 
Hermlte bicubic and the bilinear functions which are Illustrated 
In the figure and were used In the example. The final form of the 


approximate solution Is ; 


N 

v(x,t) - I q.(t) (x) . (d.4) 

1-1 


On the node points the solution is approximated by the set (q^(t)} 
for the type of hill functions which overlap as illustrated in 
Figure D-l. Solution for the set {q^(t) } is then analogous to 
solving for the set (u^(t) } on the node points using the finite 

difference technique. The Equations for the set {q^(t)} are 

» 

obtained by substituting (D^4) into (D-3). 


N 

Z 

1-1 


'0 


<^i 

3t 


*i*j + 


*1 3x3x 


+ f*.) dx 


- 0 


j - 1,2,3, N 


(D.5) 


From (D-5) a set of time differential equations is obtained which 
is solved using an implicit numerical method. 

The solution of the problem posed by the example is closely 
approximated by the erfc function in the range 0 < x < 3 if a = 6, 
and this solution was used to compare the accuracy of the finite 
element and finite differnece methods. Figures Dr2-,D-4 show the 


191 



fMUlWM* 


1 / WAI, CRIP roit.is 


Figure d-2. Error for flnltowdlfforonce ys, recfproctl of tol«l 
numbtp of grid points. 
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The error obtained In the solution by the bilinear finite 
element method Is very nearly the same as that obtained with the 
finite difference method. This was not surprising because the 
System of equations for {q^(t)} and {u.(t)} were quite Similar. 
What was surprising was that the Hermit e bicubic finite element 
produced such poor results, although this surprize was based upon 
the intuition that since this element was more difficult to 
use it should provide some reward for the difficulty. 




MICROELECTRONICS AT MISSISSIPPI STATE UNIVERSITY 
A. Introduction 

The new microelectronics facility at Mississippi State consists 
of approximately 6600 square feet of space located in the newly 
constructed (197?) Harry Charles Fleming Slmrall Electrical Engi- 
neering Building. The space and equipment represent a million 
dollars investment. Approximately 3000 square feet consists of 
class 10,000 clean rooms. Within this space, facilities exist 
for mask-making and photolithography, for chemical preparation, 
oxidation, etching, diffusion of Impurities, for metallization, 
die-bonding and lead attachment, and for the evaluation of device 
and integrated circuit parameters. In addition, the hybrid facili- 
ties allows fabrication of thick and thin film circuits and 
encapsulation. Close research ties with NASA Marshall Space 
Flight Center makes available to EE students and faculty additional 
research tools, including scanning electron microscope, Auger 
analysis, electron microprobe, etc. 

The vast field of solid state electronics has experienced 
radical changes since the invention of the transistor in 1949, and 
the Invention of the integrated circuit in 1959. Today's hand held 
calculator, for instance, contains more than 50,000 active solid 
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•tat* device* on a single substrate (called a "chip") ; electronic 
watches and appliances touch every aspect of our dally lives . 

The rapidly changing technology which has Made such electronic 
marvels possible poses a continuing challenge to an electrical 
engineering faculty not only to keep abreast, but to devise a 
curriculum which passes Such knowledge on to students, both graduate 
and undergraduate, and to motivate them In a manner adequate to 
enable them to join that technology as actively contributing 
graduate engineers. To this end, the Electrical Engineering 
Department at Mississippi State University has recently established 
an extensive laboratory facility and revised curriculum aimed at 
meeting the threefold purpose of transmitting knowledge through 
teaching , at the undergraduate and graduate levels; opening new 
frontiers through research; and extending these efforts and findings 
to others through service. 

Success in these endeavors depends on several Important 
ingredients : 

1) A dedicated faculty possessing interest and expertise in 
the various areas of physical electronics. 

2) An undergraduate curriculum devised to teach students not 
only fundamental principles of physical electronics and of solid 
state electron devices, but also of the technology which is used in 
their realization. Students should be provided the opportunity to 
achieve M hands-on" experience in device fabrication. 

3) A graduate teaching program of fering education which enables 
students to do practical design and creative research contributing 
to the state-of-the-art. 
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4) A comprehensive research program which simultaneously 
provides for the education of advanced students while contributing 
knowledge, basic and applied, to the advancement of the field. 

5) The availability of adequate laboratory facilities to 
support both teaching and research activities. 

B. New Electrical Engineering Facilities 

The Harry Charles F. Slmrall Electrical Engineering building 
was completed in January, 1977, consisting of 95,000 square feet of 
floor space at a construction cost of 4.25 million dollars. This 
building was built under the direction of the State of Mississippi 
Building Commission With funds appropriated by the State legislature. 

The general theme of the building is classrooms, computer 
terminal and other heavily populated areas on the first floor, 
department offices and Student activity areas on the second floor 
with teaching and research laboratories on the third and fourth 
floors, The features of the Harry Charles F. Simrall Electrical 
Engineering building include: 

1) 10 classrooms with 382 seats total. 

2) Computer terminal operated by computing center housing 
1004 terminal and other special purpose equipment. 

3) 254 seat Auditorium. 

4) Microelectronics laboratory complex for instruction and 
research. 

5) High voltage laboratory complex for instruction and 
research and testing. 


6) Computer design and application laboratory complex for 
Instruction and research. 

7) Dedicated laboratories for instruction and research In 
control systems, microwave systems, communication systems, and 
electronic systems. 

8) Dedicated instructional laboratories for study of digital 
devices, networks and electronics. 

9) Student activities and organization area. 

10) Faculty and Departmental office complex. 

The building now houses all programs and activities of the 
Electrical Engineering Department which have previously been housed 
in parts of Mve buildings. 

The Electrical Engineering Department includes 26 faculty 
members, some With duties outside the department, 410 undergraduate 
electrical engineering students, 42 undergraduate electronic 
engineering technology students and some 45 electrical engineering 
graduate students for a total enrollment of some 500 students. 

C. Microelectronics 

The objectives of the microelectronics program at Mississippi 
State are consistent with the University charge to provide education, 
research and service. These objectives are: 

1) To provide an outstanding facility for research dealing 
with the electronic properties of materials, electron device 
research, circuit fabrication technology, and computer aided design 
of electronic circuits. 
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2) To provide circuit technology acceeelble at several level# 
of sophistication which will support research and teaching activities 
within the department, college and University. 

3) To provide a facility available for prototype development 
of electronic apparatus by industrial concerns and to provide an 
expertise in electronic fabrication technology to assist industrial 
start-up operations. 

4) To provide operational strategies to allow utilization of 
the facilities on an economically feasible basis consistent with 
the commitment to supporting agencies. 

The SCOPE of the developing facility is detailed below; 

X. Fabrication and processing capability 
A. Silicon Integrated Circuit 

1. Thermal oxidation and Impurity diffusion 

2. Photolithography 

3. Chemical vapor deposition 

a. Polysilicon 

b. silicon dioxide 

c. other (phosphosllicate glass, silicon nitride, etc . ) 

4. Aluminum film deposition 

a. sputtering 

b. E-beam evaporation 

5. Expitaxial silicon 

6. Ion-implantation 

Hybrid Integrated Circuits 

1. Thick film screen printing and firing 
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» Chip and viva bonding 

3. Component t riming 

4. Substrate scribing 

5. Packaging and encapsulation 
C. Printed Circuits 

1. Board sensitization, exposure, and development 

2. Baord etching and drilling 

3. Board plating 

4. Board lamination 

II. Computer Aided Design 

A. Art work generation by computer 

1. Printed circuit layout from line drawing 

2. Custom art work for hybrid and monolithic Id's 

3. Art work from cell llbrairies and layout programs 

B. Art work processing 

1. Reduction 

‘ 2. Step and repeat reduction of IC master masks 

(presently thru NASA) 

3. Development of working masks 

C. Circuit Analysis 

1. PCAP 

2. Nonlinear circuit analysis program 

3. Digital logic and timing analysis 

III. Experimental Evaluation 

A. Electrical: C-V, sheet resistance, curve tracing, 
probe testing, circuit testing. 
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B. Electrochemical: lap and stain, anodisation 

C. Optical: Microscopic, interferometric, ellipsomctric 

D. Surface analysis: SEM with SIMS or Auger attachment 

(presently thru NASA) 

D. Description of Facilities 

As indicated earlier, the new $4.25 million Harry F. C, Slmrall 
Electrical Engineering Building became occupied in January of 1977. 

In the electronics area, the department owned or had custody of 
approximately $250,000 worth of equipment related to microelectronics 
prior to moving into the new facilities. Primarily the existing 
equipment consisted of that furnished by the state/ (mostly in the 
in the past 3 to 5 years), surplus equipment from NASA, much built 
in-house, and GFE items. 

The microelectronics facilities are located on the fourth floor 
of the building and consume in excess of 6,600 square feet. 
Approximately 3,000 square feet of this consist of class 10,000 
clean rooms. At a conservative building cost of $45 per square 
foot, this represents approximately a $300,000 Investment in floor 
space alone. In addition, the department either has purchased or 
is in the process of purchasing some $175,000 worth of new equipment 
to furnish the microelectronics laboratory and associated supporting 
labs. Thus, a total commitment by the University of well over a 
half million dollars is being invested in this area for purposes of 
fundamental research and teaching. In addition, equipment obtained 
through industrial donations represents well over $350,000 to date. 

203 


Through University acquisition and industrial donations, the 
department will be equipped to do significant research in silicon 
Integrated circuits, thick film and chip hybrid circuits. 

The following describes the fourth floor microelectronics 
laboratories and related support facilities. All rooms have 
compressed air, vacuum, chemical drains, and extensive venting 
for safe and convenient operation. 

A. Silicon Fabrication 

This activity is to be carried out in three different rooms, 
all class 10,000. The water supply consists of an evaporation still, 
a storage tank with ultraviolet light source, a deionizer, and 
submicron filters for delivering a minimum of two gallon/minute 
of 15-18 megohm/cm. water. . 

Alt Process Room 

This large 30 x 40 foot room has one fume hood With 
etching sink and D, I. water. Benches and cabinets and a quartz/ 
mullite tube closet are also included. There exist five vented gas 
bottle closets capable of storing ten large cylinders. Water, drains 
and 3-phase power are available for furnaces and reactor equipment. 
The basic equipment includes : 

. Therraco Spartan 6-tube furnace with class-100 load 
station and source cabinets. 

Three single tube furnaces (Lindberg, Marshall and BTU), 
all with electronic controls. 

: : . ’’ One class 100 clean bench 


. Vaponics deionized water system with additional 
reservoir lor continuous cycle filtering 
. Wafer dicer 

. Unicorp. Inc. Unipole VIII Epitaxial reactor with 
dual 8" susceptors driven by a LEPEL solid state power supply (143KVA) 
, Four point probe 
. Lap and stain facilities 
. Ellipsometer/Laser source 
, Microscopes 
. Multiple probe station 
. Capacitance-Voltage test station 
, Anodization equipment 

V Large dewars for liquid nitrogen and liquid oxygen 
For the furnaces, two tubes will be equipped to 
operate as hot-walled low pressure CVD reactors. The remaining 
tubes will be used for diffusion and oxidation. One or more tubes 
will be available for non-silicon work. 

A2. Photolithographic Room 

This room is 10 x 24 feet and is equipped 
with three in-line class 100 clean benches and two pass-through 
ovens with nitrogen ambients. All photo resist and developing 
operations are done in class 100 environment. Photo-resist 
spinning is done under temperature and humidity control. Oxide 
etching and wafer washing is done in end stations with high quality 
water. Other equipment includes: 

. Two-head spinner, plus a single-head spinner 
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. Electroglass mask aligner for 1 1/2 M wafers , plus two 
Casper aligners for 1 1/2" to 3" wafers 
. Balance table's 
, Interference microscope 
. Hot plate/stirrers, etc, 

A3. Metallzatlon and Bonding Room 

This large 30 x 20 foot class 10,000 clean room 
if for the deposition of metal films and thermal-compression 
bonding. It Is equipped with compressed air, vacuum, D. I. water 
and chemical drains. Pertinent equipment includes: 

. NRC-Varlan 6" diffusion pump vacuum station equipped 
with a Sloan sputtering system and a thickness monitor 

, Hughes pulse tip thermal compression bonder 
, Class 100 clean station for bonder 
. Varion Vacion vacuum station equipped for thermal 
evaporation of aluminum including film thickness monitor and class 
100 clean station 

, Class 100 clean benches for multiple probe station 
and vacuum system work areas 

, Work benches and cabinets 
. Large dewars for liquid nitrogen 

. A 200 KeV ion lmplanter is to be added to this room 
in early 1979. 

, Perkln-Elmer Vacion vacuum station equipped with 
dual E-beam guns for depositing alloys. 


206 


B. Hybrid Fabrication 

This activity will be carried out in one large 48 x 42 
foot room equipped with fume heads , sinks compressed air, vacuum 
and three-phase power. Ample bench tops and storage cabinet space 
is also available. Equipment Includes: 

. Thick film screen printer 
, Class 100 clean bench for printer 
. Drying oven 
. Thermco belt furnace 
. Thick film trimmer 
. Thick film trimmer 
i Multiple probe station 

, Electronic instrumentation for circuit evaluation 
. Packaging and sealing station 

C. Artwork and Mask Making 

This work will be carried out in four different rooms. 

Cl. Computer Aided Plotting Room 

This 30 x 50 foot room was initially scheduled for 
all artwork but a more comprehensive computer aided graphics 
facility is being developed. Equipment includes: 

, Gerber Scientific 1200 plotter with HP minicomputer 
tape drive (36 x 36 inch plotting surface with ink or photo exposure 
heads) 

. Tektronix graphics terminal with hard copy capability 
. Uniscope CRT terminal to the central campus UNIVAC 1108 

computer 


, Data general Eclipse S/130 minicomputer with Video 
monitors and Varlan Statoa 42 high density electrostatic printer/ 
plotter 

. Lexidata color graphics terminal 
. Digitizing table for generating check plots on 
Gerber plotter 

C2. Manual Graphics 

This will be done in the reduction room. Equipment 
Includes drafting tables and a Unitech coordlnatograph with digital 
read-out. Rubyllth masks are produced with this system. 

C3. Dark Room 

This is a small dark room with modern developing sink 
temperature controlled water source, and ample storage cabinets and 
bench top space. 

C4. Reduction Room 

This room is part of the printed circuits lab and is 
equipped with a precision Dekacon reduction camera (36 x 36 inch 
mask mounting capabilities). 

C5. Step and Repeat Room 

This is an 8 x 24 foot class 10,000 clean room and 
is part of the silicon IC complex. It is equipped with two class 
100 clean benches. One bench has a sink and hood for development 
and one bench is for a step and repeat camera. The first bench is 
equipped with compressed air, nitrogen blow gun, and a DI submicron 
filtered water. 

It is hoped that a step and repeat camera will be 

added soon. 



D. Printed Circuit Lab 

Thin laboratory la equipped with benches , a fume hood and 
sink, and two epray etchers , water dryer, KPR developer unit, 1R 
oven, dip coater, rinse unit, shea';, pilot hole punch, drill press, 
light table, etc. Almost all of this equipment is new. 

E. Other Equipment 

The microelectronics laboratory has a wide assortment of 
electronic equipment for evaluation of semiconductor devices and 
circuits, Including signal sources, meters, counters, curve tracers, 
scopes, etc. A partial listing of new equipment just purchased 
Includes : 

. 100 meg Hz Oscilloscope (Tektronics Mod 465); 

. Sampling Oscilloscope (HP Mod 182C main frame with HP 
1810A sampler); 

. Digital to Analog Converter, HP 49303A; 

. X-Y recorder, HP Mod. 7015; 

. LRC meter, HP Mod 4332A; 

. Precision multimeters, HP 3490A; 

. Interface bus, HP 59310A; 

• . 6tC . 

F. Supporting Laboratories 

FI. Electro-optic Laboratory 

This lab is equipped with a vibration free table, 
optical benches, detectors, lens, filter, optical radiation detector, 
two lasers (5 watt argon and 100 milliwatt Helium-Neon), and several 
items of standard electronic instrumentation. 
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F2. Communications Laboratory 

This facility includes extensive electronic 
instrumentation for microwave measurements and for signal analysis. 
Equipment includes wave guides, Slotted lines, microwave sources, 
UHF sources, wave analysers, antennas, etc. New equipment will 
include a spectrum analyzer (up to 1,8 G Kz) and a vector voltmeter 
for s-parameter measurements, Also included are maximal length 
pseudo-random noise generators and bit pattern generators. The 
laboratory is located in two rooms on the fourth floor and has two 
antenna locations on the roof of the building. 

F3. Instruction Laboratory 

Two large rooms and several smaller rooms are 
dedicated to teaching and project laboratories, A wide variety of 
instrumentation, breadboarding equipment, logic trainers, power 
supplies, etc. are available for instruction, training and project 
work. 
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APPENDIX F 


POST-HEAT TREATMENT EFFECTS ON DOUBLE LAYER 
METAL STRUCTURES FOR VLSI APPLICATIONS 

A< Introduction 

The Increasing demands toward greater packing densities in LSI 
and VLSI make it imperative that multilevel metallization systems 
be developed. Although several reviews have been written [Fl-7] » 

.u ' ' 

the realization of reproducible and reliable results have yet to be 
forthcoming. The most conmon problem associated with double layer 
metal has been the inability to make electrical contacts between 
metal layers through via holes etched in the dielectric to provide 
electrical communication between metals. The problem Is more acute 
when the number of vias is large (>500) and are relatively small in 
size (£0.2 mil square). Another, less Important problem is associ- 
ated with shorts between metal layers due to pinholes in the dielec- 
tric film, thin spots, or poor coverage of hillocks in the first 
level metal. An example of the use of double layer metal for cir- 
cuit application [F-8] is given in Figure F-l. 

In order to study the effects associated with double layer 
metal structures having a large number of small vias, a test pattern 
was generated consisting of a string of 560 vias. The via size was 
varied from 0.5 mil square to 0.2 mil square per string. Electri- 
cal measurements were made on the test pattern after initial 
sintering and subsequent heat treatments in order to monitor 
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Figure P-1 


An example of the use of double layer metal in the 
realization of an aluminum gate C-MOSFET structure 




contact resistance behavior and variation* in yield. Metal 
interconnect* between via* were typically 0.7 mil* wide although 
sample* were also prepared and teated having Interconnect* of 
0.4, 0.5, and 0.6 nil wide Al, with little or no variation in 
results. 


B. Preparation of feat Vehicle 

The steps involved in preparing the test vehicle consist of 
the following. The starting material was 3-8 ohm-cm, (100) oriented, 
n-type phosphorous doped silicon wafers. A field oxide was thermally 
grown in dry 0^ for 5 minutes, boiling HjO for 60 minutes and then 

dry 0 2 for an additional 15 minutes at 900°C resulting in an oxide 

■ ■ ■ 

thickness of 14 A. Prior to depositing the first layer metal, a 
cleaning step was performed, this consisted of a one-minute rinse 
in deionized water. The wafer was then dehydrated at 900°C in an 
N 2 ambient for 10 minutes. The first metal layer was d.c. sputtered 
structural grade Al alloy 6061 of 0.5ym thickness. The metal was 
patterned using conventional photolithographic techniques with 
Waycoat-31 negative photoresist. 

Prior to depositing the dielectric, scribe lines (0.7 mils 
wide) were etched around the test patterns through the thermally 
grown field oxide down to the silicon substrate. These lines were 
used as an etch end-point monitor in etching vias in the intermetal 
dielectric. 

The intermetal oxide was next deposited at 400°C using CVD of 
SiH A (4% in Ar) and 0 2 (and later P 2 0,. of approximately 3 mole 


It* 

percent) Co a thickncM of 8 A , Dielectric thlckneee was aeaaured 
with an (laser) elllpsometef, Via* vara atchad In thla dlalactric 
film with buffarad HP using standard photolithographic Making* 

Thla etching was done flrat by dipping the wafers in a stirred 
solution of etchant, using the scribe lines as reference to deter** 
mine when etching was completely through the dielectric film. In 
later processes, this etching was done ultrasonlcally In a totally 
enclosed container* 

The second layer metal wad also deposited using dc sputtering 
of the same target used in the first layer Mtal and patterned. 
Typically a 15 minute aluminum sintering process at 470*C preceded 
the first testing of the completed test vehicle. 

C. Experimental Results 

Approximately 50 wafers Were processed. Each wafer had 200-220 
test patterns with via sizes ranging from 0.2 mils square to 0.5 
mils square. Table X represents a summary of results in terms of 
percent yield obtained for each Via size as a function of fabrica- 
tion process. Process A through L are explained below* Among 
other variables it will be noted that each process Is characterized 
by an etch time associated with etching the vlas in the Intermetal 
dielectric film. The scribe lines Were used as reference. When 
all dielectric was etched from these scribe lines, this was defined 
as the 'break' time. Any additional etching beyond this break time 
is measured in seconds . 


Percentage Yield 

Via Resistance for 560 Vlas < 20 Heg Oh* 
Process Identification 

Via Size j A | 8 I 6 I 8 J E J F [ 0 I H I X i j| K] L 


0.5 nil* T 81 76 95 98 99 ] 98 94 84 94 hit 99 98 


0.4 mils I 70 56 88 94 97 98 j 92 j 82 89, 9o! 97 97 


0.3 mils 61 53 64 85 93 94 92 83 87 80 100 


0.2 mils 40 18 64 82 91 94 93 89 72 |67 


Table I-A 


Percentage Yield 

Via Resistance for 560 Vlas <10 Kll Ohms 
Process Identification a 

via size a b I I d |e|f|g|h|i[j| k| l 
0.5 mils 17 59 86 97 97 94 84 61 92 93 98 98 

0.4 mils 8 48 73 89 94 96 88 64 82 90 97 97 

0.3 mils 0 30 45 6a 89 92 68 74 85 80 100 100 

0.2 mils 0 14 |13 126 44 59 64 52 38 67l 98.| 98 

Table I-B 


Percentage Yield i 

Via Resistance f of 560 Vlas < 1 Kll Ohm 
Process Identification 

Via Size A B C D E I F G H I I j| k| L 
0.5 mils 0 59 78 85 86 82 74 56 92 92 98 98 

0.4 mils 0 46 47 66 72 72 68 54 82 90 97 97 

0.3 mils 0 28 85 36 47 55 58 56 85 80 100 100 

0.2 mils 0 8 10 15 22 27 30 30 32 67 75 98 

■ ■ ■ .-JL — -I I n i.fr I I III III.I.I III I . j ..II J 

Table I-C 

Table I Percentage yield as a function of processing with via 
size as a parameter for contact resistance of 560 vias 
less than 20 Megohm, 10 Kllohm, and 1 Kllohm for tables 
I-A, I-B, and I-C, respectively. 
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Process Definitions; 

1) Wafer* were processed as described ip lection B, 

The width of the first level metal interconnect of 0.4, 0.5, 0.6 and 
0.7 Mils corresponded to test patterns having square Vlas of 0.2, 
0,3, 0.4, and 0.5 alls, respectively. The width of the second aetal 
was 0,7 Mill, Via etch time was break plus 10 seconds. 

2) Wafers were processed a* described in section B. 

The metal interconnects on both levels were 0.7 mils wide, (This 
is true for all processes 2 through 7). Via etch time was break 
plus 10 seconds. 

3) Wafers processed as described in section B with Via 
etch time in stirred dielectric etchant extended 20 seconds beyond 
break time. 

4 thru 7) correspond to wafers processed as described in 
section B but with via etch time in stirred dielectric etchant 
extended 40, 60, 80, 100 and 120 seconds beyond break, respectively, 
(Note that the values given for processes 1) through 7) represent 
an average of five different process runs of eight wafers for each 
run.) 

8) Wafers processed as described in section B but 
having a phosphorous doped Intermetal dielectric of 3.1 mole per.- 
cent (as determined by Auger spectroscopy) and with vise etched in 
stirred dielectric etchant for 20 seconds beyond break. 

9) Wafers processed as described in 8 above except via 
etching was done ultrasonically in a closed container for 20 
seconds beyond break, 
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Yield vs» Etch Time • 




Figure F-2. (continued). Percent yield as a function of etch 
time beyond break with via size as a parameter. 
Intermetal dielectric etch accomplished with 
stirred B. 0. E. 
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10) Wafers processed as described in 9 above with the 
addition of a first level metal cleaning step through the vias 
prior to depositing the second level metal. The metal cleaning 
solution consisted of an ethylene glycol-buffered HF-l^O solution 
for the removal of A^O^ from the first level metal prior to 
deposition of second level metal [F-10] . 

11) Wafers processed as described in 10 above but having 
undergone a one-half hour post heat (sintering) treatment at 490°C. 

The results of processes 3 through 8 are displayed graphically 
in Figure F-2. The data are shown in Tables Fl-A> B, and C which 
give the percent yield for via contact resistances less than 20 
meg-ohms, 10 kil-ohms and 1 kil-ohra. It is desirable to have 
minimum contact resistance (i.e. less than 250 ohms) since as the 
number of vias increase from 600 to approximately 3000 for a VLSI 
circuit, the resistance will increase proportionately. Any pattern 
having a contact resistance greater than 20 meg-ohm is assumed to 
be an open-circuit in Table FI. 

Wafers processed as described in D in Table FI were used for 
an Auger surface analysis of the first level metal through a 0.5 mil 
via just prior to depositing the second level metal. The purpose 
of the Auger analysis was to determine what insulating matelals 
reside on the surface of the first level metal which prevented good 
low resistance ohmic contact between metal levels. The results of 
this analysis indicated the following materials were present with 
an accuracy of ± 5%: 




A1 2 0j 


64% 

C 


31% 

N 


1.8% 

Mg 


1.3% 

Si 


0.6% 

F 


0.6% 

S 


0.7% 

0 

(other than Al^Oj) 

<0.1% 


For this contact, the amounts of Mg and Si were about the same as 
the bulk A1 alloy so there did not seem to be any surface enrichment 

of MggSi which was suspected as a source of the high resistance. 

■■ . 0 ' 

When about 100A of this surface was removed by ion etching, the 
N,S,F, and C were greatly reduced, which suggested that they may 
be a source of interfacial contamination. 

The large percentage of carbon measured (31%) is believed to 
be a result of the photoresist leaving a carbon residue on the 
wafer. The composition of this photoresist could not be obtained 
from the manufacturer. 

It should be noted that the composition of the structural 
grade A1 6061 consist of the following impurities: 


Si 

0.4 - 0.8% 

Cu 

0.15 - 0.4% 

Fe 

0.7% 

Mg 

0.8 - 1.2% 

Mn 

0.15% ^ 

Cr 

0.15 - 0.35% 

Zn 

0.25% 

T1 

0.15% 


The primary reason for using this A1 alloy was for the prevention 
of hillock formation. 

Post heat treatment or sintering of a wafer exhibiting an 
initial poor yield (i.e. less than one kil-ohm for a chain of 560 
vias in series) can increase the effective yield by 500 to 700% as 
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Time of Heat Treatment (Hours) 


Figure F- 3, Percentage yield of number of chips less than 1 kil ohm 
as a function time of heat treatment at 490°C for wafers 
with undoped Intermetal dielectric. 
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Figure F-4, Average contact resistance for 560 vias of all 0.5 mil 
square chips as a function of post heat treatment time 
for wafer //2W. (Total number of chips ■ 55) 
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illustrated in Figure 3 for Wafers 10H and 2W. These wafers were 
processed as in * C * above, having only a 15 minute 470*C sintering 
before Initial testing. All addition sintering was done at 
490°C in a nitrogen ambient. 

The primary reason for this increase in yield with’ post heat 
treatment is due to the fact that the Al-to-Al contact resistance 
decreases with sintering, as demonstrated in Figure 4, Here, the 
average resistance for all test patterns having 0.5 mil square 
vias (approximately 55 of the total 212 chips) is plotted as 
a function of sintering time. As the time of sintering increases, 
the number of test patterns having a contact resistance < 1 kil-ohm 
increases, hence the yield increases. 

On an individual chip basis, the contact resistance for 560 
vias per chip behaved as shown in Figure F-5, For a given wafer 
which initially exhibits a relatively poor yield, the number of 
chips behaving in the manner shown for chips //IOC and 10W was 
approximately 51%, those behaving like chip //10R was 27% and those 
like chip // 10T was 10%, The remaining chips either exhibited low 
initial contact resistance (<300 ohms) or were open-circuited 
(>20 meg-ohms). Those chips initially exhibiting a low contact 
resistance generally had their resistance lowered with increasing 
heat treatment, how4ver a few chips showed a gradual Increase in 
resistance with Increased heat treatment time as shown in Figure^F-6. 
It is believed that this increase in resistance is due primarily to 
the consumption of A1 by SiO, to lessen the A1 thickness, Since 
at 490°C A1 reacts chemically with Si0 2 to lessen the A1 thickness. 
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Since at; 490*C A1 reacts chemically vith SiOg to form and 

0 

free Si,Al is consumed at a rate of 210A/hr. [F-15] 

The most predominant trend found in sintering of double level 
metal test patterns is demonstrated in Figure F-7. AS inferred 
from this plot of contact resistance as a function of post-heat 
treatment for these two wafers* the average contact resistance 
decreases drastically for the first few hours of sintering* 


Because of this feature* the realization <vf multiple level metal 
LSI having a large number of vias (>2000) is a distinct possibility. 

One major problem involved in post heat treatment of double 
layer metal structures is that the dielectric (C.V.D. - SiO^) has 
a tendence to crack or "craze" at elevated temperatures (>500°C). 
This cracking, as illustrated in Figure F-8, is due to excess 
stress on the Insulating layer [F.9]. One method to lessen, if 
not-alleviate this problem, is to use a phosphorous doped dielec- 
tric of 3 to 5 mole percent. This increases the effective tempera- 
ture for which sintering can be done before cracking is observed [F.9]. 

Using this phosphorous doped dielectric (also called phos- 
phosllicate glass or PSG) resulted in a much higher yield as shown 


in Table 1 and also a considerable less contact resistance as 
shown in Figure F-9. Here the 1 average contact resistance as a 
function of via size is plotted for a test pattern with a doped 


dielectric layer for no heat treatment (curve BO and with 30 
minutes post heat treatment at 490°C (curve A) . These results are 
compared to the best results obtained for the undoped case (curve'C) , 
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Contact Resistance of 560 Vlas (Meg-ohms) 



Figure F~5, Contact resistance of 560 vias as a function of heat 
treatment time for several different chips of wafer 
#10H at a temperature of 490°C. 

225 




Contact Resistance of 560 Vlas (Ohas) 



0 1 2 . 3 4 5 6 7 

Time of Heat Treatment (Hours) 


Figure Fie. Contact resistance of 560 vlas as a function of heat 
treatment time for wafer //10H showing slight increase 
In resistance with time of heat treatment for -6% of 
the chips on the wafer. 
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even after many hours of heat treatment! The table in the 
figure Indicates resistance and percent yield for each case 
presented. 


A resistance map of the wafer shown in curve B Is given below; 
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D. Analysis and Conclusions 

The following comments and conclusions may be derived from 
this investigation. 

1) The variation in width of the metal Interconnects between 
Vlas had little influence on yield, as expected. It may have had 
a slight effect on total resistance, expeclally for the 0.4 mil 
compared to 0.7 mil wide interconnect, but this was not detectable. 
The major difficulty came in mask alignment of the smaller Interconnects. 
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Average Resistance of 560 Via Contacts" (P his ) 


700 


Figure 



F-9. Average contact resistance vs. via size for wafer 014-1 

having phosphorous doped Intermetal dielectric before post 
heat treatment (curve B) and after 30 minutes at 490°C (curve A). 
For comparison with a wafer having an undoped dielectric, curve C 
represents such a wafer after undergone 8 hours of post heat 
treatment at 490 # C. (Note contact resistance for 0.5 mil via 
of wafer in curve C prior to post heat treatment was 90 kilohms.) 



2) The extension of the via etch time of the dielectric film 
for both the dipped (stirred) etching and the ultrasonic etching 
indicated a substantial Improvement in yield with only a slight 
undercutting for extended times below 90 seconds. (As expected, 
the doped dielectric etched faster than the undoped.) For etch 
times beyond break in excess of 110 seconds; the yield decreased. 
This is believed to be attributed to an Increase In contact 
resistance- and thus a decrease in yield - due to a chemical 
reaction between the dielectric etchant and the alloys of the 
sputtered aluminum to form an insulating layer at the surface of 
the first level metal prior to depositing the second layer. 

As Indicated in Figure F-2, the optimum etch time is dependent 
on the size via being etched. 

3) The addition of phosphorous doping to the oxide eliminated 
the possiblity of intermetal dielectric cracking for post, heat 
treatment temperatures below 500°C. A considerable increase in 
yield was also realized, having total average contact resistances 
slightly over 200 ohms for 0.5 mil vias and 250 ohms for 0.2 mil 

as per Figure F-9. (A slow pull of the wafer from the furnace was 
also practiced to eliminate cracking. ) 

4) The use of an ultrasonic etching bath for etching vias 
gives a much better probability of etching small vias (<0.2 rails). 
This etching was accomplished in a completely enclosed container. 

5) The chemical cleaning of the first level metal through 

the vias prior to depositing the second level metal also increased 
yield, The first layer A1 develops a thin oxide coating as soon as 
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it is exposed to the atmosphere and this oxide skin creates a 
high resistance upon depositing the second layer. This coupled 
with the difficulty of removing all of the dielectric film from 
the via and the formation of oxides of the alloys contained within 
the Al all tend toward low yield, high resistance contacts. There 
are actually two approaches for solving this problem, back sputtering 
and chemical etching. Back sputtering has the advantage that the 
oxide thickness can in principle be reduced to zero if the subse- 
quent metallization is applied without exposure to air. The 
disadvantages include: possibility of radiation induced MOS damage; 

potential for contamination of the sputtering system; redepositlon 
of previously sputtered materials especially from the substrate 
table; etc. The chemical etch procedure does not have these 
disadvantages, but the oxide thickness can never be reduced to 
zero. [F. 141 

The use of a 1:1:1 ethylene glycol: buffered -HF : 1^0 etch 

O 

removes all but 30-50 A of oxide [F. 10] without attacking the Al. 

(The second layer sputtered Al can easily penetrate this oxide 
thickness to form good ohmic contacts.) This etch was done just 
prior to second level metallization. 

6) The application of post heat treatments can drastically 
Increase yield by lowering total via contact resistance. There 
are Several possible reactions which could occur to cause this 
temperature effect. [F.ll] Aluminum could react with aluminum 
suboxides and thus break up the continuous layer of dielectric. 
Alternately, the recrystallization of aluminum in both of the 
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level could mechanically disrupt the thin native oxide layer, 
therefore forming aluminum-to-aluminum contact. Ara yet, there 
is no evidence for this speculation. 

Figures F-4 and F-5 Indicate that In some cases the contact 
resistance Increases with post heat treatment time prior to the 
"diffusion" through this 'interfacial contamination' by the two 
layers of aluminum. If this contamination were SiO^, then It 
could be said that there exists excess ionic silicon in the oxide 
which becomes tled-up upon the application of post heat treatments 
[F-l]]. Since the contamination is probably Al^O^, one might say 
there exist excess ionized A1 in the oxide which reacts with the 
oxidizing species during the sintering process, thus decreasing 
the total ionic charge in the oxide layer causing the net resistance 
to increase. 

7) Testing the via test pattern can be tedious. For high 
resistance vlas the measured resistance was seen to be both light 
sensitive and current sensitive. As a result, all testing was 
accomplished in the dark. Ideally, very low testing potentials 
are applied In order not to break down a barrier layer thus 
turning a defective via into a good one* Also, theoretically, 
vlas should not be tested in series, since the entire applied 
voltage will be dropped across a defective via, probably causing 
it to break down and appear good. However, in our testing, all 
Vlas were measured at once in series using a digital voltmeter. 

8) Recently, a new planar multilevel interconnection tech- 
nology was introduced [F. 13] using polyimide films. Although 


I ' 

thli procedure appears to be very promising for future applications 
If taken at face value. It Is felt there Is still many facets of 
M maglc M associated with It. 
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